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ABSTRACT 

Two models for the correlation of ship position reports obtained 
and relayed by ocean surveillance satellites are described. The models 
are intended to process position reports of ships detected by satellites 
and correlate this data to provide intelligence concerning the position, 
course, and speed of ships traversing the ocean surface. A probabilis- 
tic model assigns position report correlations by reliance on probabil- 
ities assigned to computed ship speed and to changes of course and 
speed. <A reliability factor is determined for each track. The deter- 
ministic model avoids explicit use of probability theory and correlates 
position reports on the basis of proximity of the position report to a 
position predicted by track extension. Evaluation of the models by means 
of computer simulations is described. Results of the evaluation are 
discussed and recommendations are made concerning areas for further 


investigation. The computer programs are included in the Appendicies. 
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1. Introduction 

Since the first sputnik was successfully projected into orbit in 
1957, the possibility of using a satellite for reconnaissance and 
surveillance of the oceans has been of increasing interest to military 
men. Studies have been conducted to analyze the usefulness and feasibil- 
ity of a surveillance satellite system and several reports have proposed 
system requirements, A report by the Pacific Missile Range described 
some of the problems to be encountered in processing data received from 
a surveillance satellite. (2] 

During the spring and summer of 1964, analysts at the Advanced 
Programs Division of the Pacific Missile Range conducted a computer 
simulation study employing a surveillance satellite to detect and report 
the positions of ships at sea. The purpose of the study was to determine, 
for a given satellite sensor capability and orbit, the times and positions 
of ships at sea eligible for detection by the satellite. The ship data 
used in the study was extracted from actual weather reports. (The studies 
are reported in detail in Pacific Missile Range Technical Note No. 3285- 
570, 1964.) The Pacific Missile Range program permitted the user to 
describe the orbit of a satellite in altitude and inclination to the 
earth's equator, It also permits the analyst to specify several different 
sensor capabilities. Ship's call signs, reported position, time, and 
cloud cover were extracted from the weather reports, and from this, tracks 
were projected. The simulation program then determined the times and 
positions of the ship at which the satellite would be in a position to 
detect the ship with each of the sensors described by the user. The 
simulation program also determined whether this conjunction of ship and 
satellite occurred during daylight or darkness. A Monte Carlo process 


based on reported cloud cover was used to determine whether the ship 
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was obscured by clouds at the time of conjunction. 

The object of this thesis is to provide the next logical step in 
the analysis of a surveillance satellite system, i.e., to devise a 
correlation logic to determine to what extent the sets of conjunctions 
provided by the Pacific Missile Range simulation can be processed to 
form tracks which correspond with the actual track as determined from 
the ship's weather reports. It will then be possible to determine the 
coverage frequency required to provide an acceptable tracking reliabil- 
ity, hence the number of satellites required. The authors have there- 
fore addressed themselves to the task of devising correlation logic to 
provide programs which will facilitate further study of the feasibility 
and cost of tracking ships at sea by use of data collected by surveillance 


satellites. 





2. General Analysis 

The formulation of track correlation logic which will facilitate 
testing the feasibility of satellite surveillance of ocean shipping, 
has been undertaken via two different approaches. The two procedures 
considered are described in detail in the following two chapters. In 
the development of both procedures it is assumed that the input data 
to the tracking program will consist only of the time of conjunction 
and the geographical position of the contact at the time of conjunction. 
Thus the analysis is not dependent upon the development of a high reso- 
lution sensor which might provide information about the ship's saaeee 
by virtue of its orientation to the orbital path, or other identification 
clues, such as length or the number of stacks observed. 

This study has been limited to a consideration of methods by which 
the data described above might be correlated in a computer program to 
derive tracks of ships at sea. No attempt has been made to determine 
optimal orbits or combinations thereof, 

It is evident that ocean traffic intelligence gathered from other 
sources would be of significant value in decreasing the uncertainty 
involved in surveillance by satellites. Such other sources include 
sighting reports by aircraft and friendly shipping, reports of ship 
arrivals and departures gleaned from commercial sources or news media, 
and voluntary reporting of own position by ships at sea. If one half 
of the ships in a certain area can be identified by means of external 
intelligence, the number of possible correlations to be determined may 
be reduced to one quarter of what would otherwise be required. This 
report will not concern itself with the methodology or effect of such 
an input to the tracking program, for if tracking can be shown to be 


feasible using information gathered solely by the surveillance satellite, 
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it will clearly be enhanced by inclusion of externally gathered 
information, 
2el The Correlation Problem 

The correlation problem can be thought of as consisting of two 
parts. The first and most difficult part involves track initiation, 
i.e., making a decision as to which pair of positions reported on 
consecutive orbits form the most probable true correlation. Once this 
decision has been made, the track can be considered as true and the 
course and speed of this designated track then becomes available for 
the second part of the problem, the track extension. 

Track Initiation. The speed with which ships sail the surface of 
the oceans is less than some generally known maximum speed. Few ships 
today can sustain a speed in excess of 40 knots for more than a few 
hours. An analysis of merchant ship cruising speeds reveals that their 
distribution is approximately normal with a mean of 14.0 knots and a 
standard deviation of 2.5 knots. Only 2.0% of the cruising speeds ex- 
ceeded 19 knots while 78.5% were in the range 11.5-19.0 knots, The 
inclusion of military ships in the analysis presumably would bias the 
results toward somewhat higher speeds but even fast naval ships seldom 
cruise at speeds in excess of 20 knots. Thus it may be reasoned that 
a correlation between two position reports which indicates a speed of 
16 knots is much more likely to be a true correlation than one which 
indicates a speed of 10 knots or 25 knots, The existence of this speed 
limitation determines a selection rule for track initiation which per= 
mits the elimination of all possible correlations between position 
reports which would require a speed in excess of some specified maximum 
speed. In many cases this selection rule alone will determine the cor- 
relation. i.e., if for a point from one sweep there is only one 
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corresponding point from the next sweep which is within a circle whose 
radius iS equal to the maximum speed attainable multiplied by the time 
between the position reports, then there is a high probability that 
these two points can be correlated as a true track, It is evident that 
any tracking program must not deny correlations showing very high or 
very low speeds, but much more confidence can be placed on a correlation 
which indicates a speed in the range of 12-18 knots. 

The difficulty of correlation increases when a number of possible 
correlations occur within this maximum distance circle. The number of 
tracks to be considered as possible for position reports from two orbits 
is equal to the number of points from the second orbit which lie within 
the maximum distance circles drawn about the points from the first orbit. 
This situation is illustrated in Figure 1. In aregwof high shipping 
density it is to be expected that the picture will be much more compli- 


cated than that illustrated, 





Note: The notation p,,, refers to report number n on orbit number r. 
Illustration of Possible Tracks for a Given Situation 


Figure 1 
5 





When the circles of maximum attainable distance are drawn in 
accordance with the selection rule for track initiation it can be 
seen that while Pj, may link with all four points obtained on orbit 


two, the same is not true for p and Py, » Thus it can be 


12° P33” 


assumed that since Py, can only correlate with Po, 


then cannot 
’ Pai 


be logically considered as a possible correlation for Pay - A similar 
reasoning applies to and eee is removed from considera- 
§ app Pi3 ’ Po) Poy, era 


tion with then can be linked with - the result of ly- 
Pi, Pi > Po3 apply: 


b 
ing the selection rule on the position reports illustrated in Figure 1] 
therefore reduces the eight possible tracks illustrated to the following 


firm correlations: 


‘ee 28 
ake) A 
i) Weak 
‘ii 22 


This selection rule cannot stand alone however, Since it is 
heavily dependent upon the assumption of a sensor which is 'perfect! 
in the sense of a radar with a blip-sScan ratio of one, i.e., if a 
ship is within the sweep path of the sensor it will be detected by the 
sensor with probability one. It cannot be strictly applied when the 
circle of maximum attainable distance drawn about a point encloses a 
seaport since there iS8 a possibility that the ship did enter port and 
was not detectable by the satellite. Similarly, any ship detected within 
a circle of maximum attainable distance drawn about a seaport may be a 
ship just starting her voyage from that port. The rule also ignores the 
possibility of a temporarily surfaced submarine or occasional data points 


resulting from contacts with aircraft or heavy clouds. 





Track Extension. After tracks have been initiated, a rule of 
consistency utilizing both course and speed can be applied to determine 
the extension of tracks for position reports from consecutive orbits. 
This rule is one that states the extension of initiated tracks to new 
position reports can be done on the basis of course and speed consis- 
tency. Ships at sea are generally steaming from one port to another 
port and have selected a direct course. This implies that most course 
changes in general are executed to circumnavigate points of land or are 
occaSioned by circumstances such as avoidance of severe weather, mid- 
course changes in destination, or navigational course adjustments 
necessitated by the effects of wind and sea. Ship cruising speeds are 
also usually constant and are determined by economical considerations, 
with slight variations to permit favorable entry conditions at the 
destination. 

There are exceptions to this rule of consistency. Fishing vessels 
off the Grand Banks follow very erratic maneuvering patterns. Navy ships 
executing training maneuvers or on certain patrol missions do not follow 
consistent course patterns. Since these exceptions apply to only a 
small fraction of world ocean traffic, a sequence of data points, i.e., 
position reports, which indicates a consistent course and speed is 
probably an actual ship's track. 

There are certain geographical areas where the antithesis of con- 
sistency is true. For example, ships approaching the Straits of Malacca 
from the north-east will, with a very high degree of certainty, change 
course radically to the north-west at the mouth of the strait. A highly 
sophisticated correlation program will be cognizant of the conditions 
and locations where these exceptions occur. However, even lacking this 


sophistication, much of the confusion resulting from circumnavigation 
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course changes will be resolvable by lending more credence to speed 
consistency than to course consistency in these areas. 
22 Application 

Each of the selection rules discussed have certain application 
limitations. One must conclude that an optimal correlation procedure 
would use them both to a certain degree, or in certain cases, but be 
able to distinguish when the selection rules break down, It is also the 
case that the longer the sequence of positions to be considered by the 
correlation logic, i.e., the smaller the interval between orbits, the 
more confidence one may place in the results. Certainly, there can be 
no less than a sequence of three points for consideration in order to 
initiate and extend tracks. 

The two correlation programs considered by the authors are described 
in detail in the next two sections. These programs employ both of the 
selection rules to some extent, but each emphasizes a different aspect. 
The first program, the deterministic correlation model, examines several 
sets of position reports at one time and makes maximum uSe of the selec- 
tion rule for initiating tracks, i.e., the elimination of consideration 
for correlation of two position reports from consecutive orbits which 
are separated by a distance greater than the distance attainable by a 
ship cruising at the arbitrarily stated maximum Speed throughout the 
period between orbits. The second program, the probabilistic correlation 
model, emphasizes the selection rule for track extension, i.e., the 
greater reliance on consistent course and Speed to select the most 
probable correlation. After considering each new data set, this program 
presents a list of all possible tracks together with a credibility factor 
which reflects the consistency of course and speed indicated by each, 


The probabilistic model will be described in detail in section four, 





Section three contains a discussion of the logic, step-by-step proce- 
dure for correlation, and details the methods and results of the 


simulation used to evaluate the deterministic correlation model. 





3. The Deterministic Correlation Model 

The deterministic model attempts to correlate reports of the 
position of ships at sea as obtained and relayed by surveillance satel- 
lites. Explicit use of probability theory is avoided. No reference is 
made to information such as normal cruising speed, or the existence of 
the more or less well defined shipping lanes between seaports. External- 
ly gathered intelligence information concerning shipping is not used by 
this model in its present level of development. The model is designed 
to permit track initiation and track extension by accumulating position 
reports from several orbits and correlating them concurrently. As it 
now stands, this model is a mid-ocean model. It has no logic for cor- 
relating position reports in proximity of seaports. 

To facilitate formulation of this model, certain assumptions have 
been made about the capabilities and characteristics of surveillance 
satellites and ocean shipping. 

These assumptions are: 

1) On every orbit the satellite borne sensor will report the position 
of all ships on the ocean that are within its range of detection. 

2) All positions reported will be positions of ocean shipping. 
Clouds, aircraft, land and objects on or over land will not be reported. 

3) Two positions reported at times differing by dt cannot corre- 
late if they are separated by a distance greater than the product of dt 


and S$ 


ees where a is an arbitrarily stated maximum speed attainable by 


ocean shipping. 

4) A ship reported on orbit i will also be reported on orbit i+ 1. 
That is, the ship cannot enter port and thus avoid detection. This dis- 
regards the possibility of submarines being detected on orbit i and sub- 
merging before orbit i +1. 
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3.1 The Correlation Logic 


Let P be t t of ti orate 
bag e the set of n positions, Psa , Pp Dp 
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ot aes i,p 

reported on orbit i. Let Q; j be the set containing all of the positions 
9 

off set Pitlym , which are within a distance d of Pij where d is the 

maximum distance a ship can travel in the time between orbits i and i+1 

at speed 5. , and Pi j is an element of P; , . Define ak the number 


9 


of elements in the set Q. ee Then 
9 


ion” Wi, iO aa a al a Psa} 


Qs = }Pia2 x | Pita © Pid m and is within distance d of Pa, 


The number of elements in Qj is greater than or equal to one, i.e., 
1£n,=2, since under the assumptions stated above, all ships reported 
at orbit i will also be reported at orbit i + 1 and none will have travelled 
a distance greater than d, the maximum attainable distance at apeed Ox ° 

The step=-by=step procedure for the concurrent correlation of ship 
position reports in the deterministic model is: 

1) Construct the sets Qs for all Desig 

2) Determine which, if any, of the sj have only one element, 

3) If any position which is the only element of a one element set 
is also included in any other set Uj , remove it from those other sets. 

4) Repeat steps 2 and 3 until no more changes can be made. 

5) For each Qj containing only one element a firm correlation 
exists between a position of orbit i and a position from orbit i+. 
Confirmation of this correlation is sought by projecting the track to the 
next (or prior) orbit. If the track is confirmed by the presence of a 
ship in the area of the projected position, the correlation is designated 


as a track, 
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6) Attempt to trace each designated track from the time of its 
designation back to orbit one and forward to the last orbit reported. 

7) Attempt to correlate tracks for the sets ar containing more 
than one element, 

8) Attempt to correlate all remaining uncorrelated position reports. 

Step 6 above is accomplished by the following methodology: 

a) Using the course and speed indicated for the designated 
track, determine the next projected ship position by projecting the track 
ahead (or back) to the time of the next (or previous) orbit. 

b) Define an area about this projected position and determine 
whether an undesignated position report from the next (or previous) 
orbit is located within this area. The size of the area defined can be 
varied as a function of the expected inaccuracies of reported positions. 

ec) If only one ship position report is in the defined area, 
and if that position report is not simultaneously in the area defined 
about another projected position, it is assigned as an extension of the 
designated track, This position report is then removed from all other 
sets, 8 , which may contain it. 

d) Repeat (a) through (c) until each designated track is 
extended as far as possible. 

e) Return to Step 2 until no further designations can be made 
on the basis of the deletion of elements from the sets Si j ; 

If more than one ship position report, or if no ship position 
report, is in the defined area, the projected track cannot be extended 
by the methodology for Step 6. 

At the completion of Step 6, the methodology for Step 7 is: 


a) For each Q; j consider the relation between Pi ij and 
bd 9 
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each element of WO as a tentative track. 

b) Determine the projected positions for each of these tracks 
by projecting each tentative track to the next (or previous) orbit. 

c) Define areas about each of these projected positions. If 
only one of the areas contains a position report, this position report 
is considered as defining a designated track. 

d) For all tracks so designated, return to Step 6. 

Those Uj which still contain more than one element, as a result 
of the areas constructed about the projected positions containing other 
than one ship, are analyzed to determine which ship is closest to the 
projected position. That position report which is closest is considered 
to be the most likely correlation and is designated. The track is then 
extended as in Step 6. 

In summation, the correlation logic for the deterministic model is 
based upon the previously stated assumptions and upon the following 
concepts; 

1) Any set 8% 3 which contains only one element determines a unique 
correlation with certainty of being correct. 

2) Tracks which are tentatively extended to a constructed area 
which is occupied by more than one ship, can be correlated on the im- 
plicit probability that the ship nearest to the center of the area is 
the most likely true correlation. 

3) Any track projecting to an area which contains no position re- 
ports does so as the result of a ship course or Speed change. Based upon 
the implicit probability that a small course change or a small speed change 
is more probable than a larger course or speed change, the undesignated 


position report nearest to the center of the constructed area is considered 
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to be the most likely true correlation. 

Although some of the correlated positions may have been correlated 
erroneously due to changes in the course or speed of the ships, the 
assumption that each ship is reported on each orbit assures that event- 
ually all positions will be correlated. 

3.2 Applications of the Model 
A simple situation illustrating one of the basic principles of this 


logic is diagrammed in Figure 2. 


Qir 





Illustration of One of the Basic Principles 
of the Deterministic Correlation Model Logic 


Figure 2 


Four position reports from orbit one and four position reports 
from orbit two are plotted and the Q. j have been constructed about 
9 


each P The true tracks are drawn in solid lines with other poss- 


ee 
ible tracks indicated by dashed lines. The position report p) j can 
b 
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be correlated only with members of the associated Q . To permit other 
b 
correlations would violate the assumption concerning the maximum speed 


attainable. Due to the assumption that the position of each ship reported 
at orbit one will also be reported at orbit two, we are assured that each 
of the eight position reports plotted can be correlated. Therefore, the 


next step to be taken after the sets Q, , have been constructed, is to 
b 


determine which, if any, of the Q have only one element. From Figure 2 


i) 


it can be observed that 21 3 and oy 


mines a correlation between Py 3 and P. I and a correlation between Py ; 
9 9 9 


have only one element. This deter- 


and Oe 
It is clear that a specific ship cannot occupy more than one position 
at a given time. Therefore Poa and Poy can be removed from set 2p 


and Po 4 can be removed from set Q Further examination now shows 
> 


ee 
that % 2 has only one element and a correlation can be established be- 
tween Bing and Deas After deleting Po 3 from set ae a correlation 
is determined between Pia and Po 2° 

Although Pia had four tentative tracks, a unique correlation of 
Pal and = has been made. This simple, easily resolvable situation 
is not likely to be encountered frequently, but it does Serve as a basis 
and as a first step toward resolving more complicated situations. 

Figure 3 illustrates a situation requiring theapplication of all 
of the steps outlined in section 3.1 above. Four position reports from 
each of four orbits are plotted. Actual tracks are indicated by solid 
lines while tentative tracks and projections thereof are indicated by 
dashed lines. The circles indicate the areas constructed around pro- 
jected positions. The circles defining the sets Q; j are omitted for 


9 


the sake of clarity. 
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The Q) j are constituted as follows: 
b] 


1.17 Par? Poa’ Poy 
ee ae > Py, 
“yee "ome ©2035 22) 
M147 Par? Poa Poy 


The tentative tracks contained within “5 ; have been omitted from the 
figure for clarity. 

Since none of the Q contain single elements, the tracks are 
initiated by projecting .. of the tentative tracks and checking for 
position reports in the constructed areas. In this case only one of 
the tentative tracks projected from each of the P, results in an 
occupied area. Therefore all four tracks are initiated correctly. When 
track one and track three are projected from orbit three to orbit four, 

a conflict situation occurs. Both projections find an area occupied by 
two ships. This situation is resolved by comparing the distances of the 
tentative projections from the center of the area projected by the course 
and speed of the designated tracks, i.e., track one is projected to orbit 
four and finds Pad and P, 3 both in the area. However, Prd is nearest 
to the center of the area and is the favored correlation. 

Tracks two and four, when projected from orbit three to orbit four, 
find only one ship in their respective projected areas, but it is the 
same ship, a A comparison of the distance from P, 2 to the center 
of each of the projected areas shows track two to be favored and leaves 
track four with an empty area. The program must then look to the other 
members of 93 4° Since, due to the deletion of elements from 83 4 as they 


were asSigned designations by other tracks, Pry is the only remaining 
b 


element in Q, 4 it is correlated with track four. 
bd 
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Diagram of Track Initiation and Extension 


Figure 3 





These two illustrations by no means point out all of the con- 
plications that may arise in track correlation. They are intended 
merely to demonstrate some of the fundamental points of the logic. 

3.3 Model Evaluation 

This correlation model was programmed for evaluation on the CDC 
1604 computer located at the U.S. Naval Postgraduate School computer 
facility. The program was written in the FORTRAN 63 computer language. 
The FORTRAN 63 program is listed in Appendix I. 

The program was designed to have the capability to receive a wide 
variety of inputs. Data representing specially designed shipping 
Situations can be entered into the program. Alternatively, situations 
can be generated by the computer program by use of random number genera- 
tors. The program is written to correlate positions gathered from 
several orbits concurrently with none of the tracks designated prior to 
the initial orbit. It is considered that ability to resolve such a 
situation is more demanding than a situation where the program is merely 
asked to extend tracks previously designated and correlate them with new 
position reports. 

Parameters that can be varied to permit full analysis of the corre- 
lation program capabilities include: 

1) The number of ships to be tracked. 

2) The number of orbits to be processed concurrently. 

3) The time between orbits. For correlation of positions in a 
concentrated area, the speed of the satellite passage is so large rela- 
tive to the speed of the ships being tracked and the area concerned, 
that it does not invalidate the evaluation to consider that all of 
the positions reported during one orbit over the target area have been 


taken simultaneously. 
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4) The size of the area containing the initial positions of all 
ships. This parameter determines the shipping density. 

5) The distribution of initial courses. 

6) The distribution of initial speeds. 

7) The distribution of initial positions within the playing area. 

8) The distribution of the number of course and speed changes to 
be made. 

9) The distribution of the magnitudes of course and speed changes 
to be made. 

10) The size of the areas to be constructed about the projected 
positions. If it is required that the area size be incremented to look 
for ships that have changed course or speed, the program specifies that 
both dimensions of the area will be doubled, | 

Provisions are also included that permit bypassing the routine 
which looks for the sets Qs; which have only one element. When this 
provision is bypassed, the tracks are initiated by projecting all 
tentative tracks ahead to the next orbit and checking for area occupancy. 
Decisions are then made on the basis of minimizing the distance between 
the center of projected areas and the correlated positions. 

The first evaluation of the model involved the correlation of ten 
position reports received per orbit for five orbits. The time between 
orbits was chosen as one hour. The initial positions of the ten ships 
were uniformly distributed in a 200 mile square. Initial speeds were 
uniformly distributed in the interval 5-35 knots. Initial courses were 
uniformly distributed on the interval 0-360 degrees subject to the con- 
straint that no ship was assigned a course which would prohibit that ship 
from being included in the ij associated with a position of some other 


ship. No course or speed changes were made. 
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This situation was correlated once using the entire logic and 
once omitting the check for Qi j which contain only one element. All 
positions were correctly correlated, designated and tracked for all 
orbits. There was no detectable difference between the evaluation with 
and without the check for single element Qs 3 P 

The above situation was modified by changing the area containing 
the initial positions from a 200 mile square to a 100 mile square, and 
the evaluation was repeated with identical results, 

Using the same initial positions, initial courses and speeds, and 
the same area, the situation was modified to provide large course and 
speed changes at the time of orbit three for three ships. Additionally, 
one ship was provided with a new course and a new speed on every orbit. 
The model was again successful in correctly correlating, designating, 
and tracking all ships thru all orbits. 

The model was next evaluated using randomly generated data under 
the following constraints; 

1) Ten ships were to be tracked. 

2) Five orbits were to be processed concurrently. 

3) The initial positions of all ships was uniformly distributed in 
an area by use of a random number generator. The size of the area was 
varied as a parameter. 

4) Initial courses were generated from a uniform distribution on 
the interval 0-360 degrees. 

5) Initial speeds were generated from a normal distribution with a 
mean of 15 knots and a standard deviation of 5 knots, subject to the con- 
straint that the minimum speed was 5 knots and the maximum speed was 35 


knots. 
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6) For each of the last three time intervals the number of ships 
to change course and the number of ships to change speed were generated 
from a uniform distribution on the interval 2-8. 

7) The determination of which ship was to make each change was made 
by generation of a random integer on the interval 1-10. 

8) The magnitude of the course change was generated from a normal 
distribution with mean zero and with a standard deviation of 0.50 radians. 

9) The magnitude of the speed change was generated from a normal 
distribution with mean zero and with a standard deviation of five knots 
subject to the constraint that the resulting speed must lie on the inter- 
val 5-35 knots. 

10) Given a designated or tentative track at orbit i, the initial 
area to be searched for the extension of this track is centered on its 
projection from orbit i to orbit i+1l. The area extends + 2T miles 
radially and = 0.09 radians in azimuth, where T is the time interval 
between orbit i and orbit i+ 1. The radial dimension of the area is a 
function of T, while the angular dimension is a function of the distance 
from the area center to the track position at orbit i, i.e., a function 


of track speed and T. Figure 4 illustrates this area definition: 


yT 
& a al 
hr qdinv 
+a -----F5 [Re 
Track At TRACK AT J 
sl ea orbit L a iecte at 
wos te 
ov 


Definition of Area Constructed Around Projected Positions 


Figure £4 
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It is to be noted that the number of course and speed changes 
permitted for this test significantly exceed the number to be expected 
under any but the most peculiar situations to be encountered at sea, 
Further, the standard deviation stated for the magnitude of course and 
speed changes is considerably larger than that which would be expected 
under normal conditions at sea, These larger values are specified for 
this evaluation to provide a more severe test of the capabilities of 
this model. 

Figure 5 summarizes the test simulations, the simulation results, 
and indicates the parameters which were varied. A sample simulation 
and the results of the correlation are contained in Appendix II. While 
the availability of computer time prohibited conducting sufficient 
simulations to permit a rigorous statistical analysis of simulation 
results, the number of simulations conducted did indicate some general 
relationships between the various parameters and the results. A com- 
parison of columns five and six of Figure 5 indicates that a major por- 
tion of the incorrect track extensions involved the extension to the 
position reports obtained on orbit one. This is attributable to the 
much greater shipping congestion at that time due to the initial posi- 
tioning constraints. A major drop in program performance is indicated 
for the evaluations with an orbit interval of 8.0 hours and an initial 
area of 50 miles square or of 100 miles square. This drop is explained 
by noting that for each of the 80 simulations conducted under these con- 
ditions, each of the Qf constructed about the position reports from 
orbit one contained all of the positions reported on orbit two. A modi- 
fication to preclude any course or speed change until three orbits had 


been completed was made for 20 simulations at 8.0 hour orbit intervals. 
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This required all ships to maintain a steady course and speed for 

16 hours. Under these conditions the effectiveness increased from 49% 
to 72% as indicated in Figure 5. In all of the simulated situations, 
many of the ships that were not correctly tracked through all orbits, 


were correctly tracked through four orbits. 
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Column l Column 2 Column 3 Column 4 Column 5 Column 6 


SIZE OF SQ. AVERAGE NO. PERCENTAGE PERCENTAGE 
CONTAINING OF COURSE NUMBER ORBIT OF SHIPS OF SHIPS 
INITIAL AND SPEED OF INTERVAL CORRECTLY CORRECTLY 
POSITIONS OF CHANGES PER SIMULATIONS (HOURS) TRACKED TRACKED 
ALL SHIPS SIMULATION THRU ALL THRU LAST 
ORBITS FOUR ORBITS 
100.0 26.9 42 220 94.0 96.2 
50.0 eins id, LA 82.0 90.5 
50.0 16.3 22 2.0 85.4 93.7 
50.0 23.8 4,0 8.0 28.3 86.5 
100.0 23.0 4,0 8.0 49.0 90.7 
*100.0 DOD 20 8.0 72.0 93.5 
500.0 29.8 20 8.0 81.5 90.5 


* For these simulations, all course and speed changes were made for 
the last two orbits. 


Summary of Test Simulations for Deterministic Model 


Figure 5 
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3.4 Conclusions 

The following conclusions are permitted by the results of the model 
evaluations discussed in section 3.3 above; 

1) Under the stated assumptions, ships that maintain a steady 
course and speed will be correctly correlated, designated, and tracked 
in almost every case. 

2) For a shipping density of 10 ships per 50 mile square or less 
and for an orbit interval of 2 hours or less, the model performs track 
correlation with a high degree of reliability and precision, even when 
a large number of course and speed changes are permitted. 

3) For an orbit interval of 8.0 hours the model correlates tracks 
with a high degree of reliability and precision after the ships have 
departed from the congested area imposed by the initial problem condi- 


tions. 
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4. The Probabilistic Model 

The probabilistic model attempts to correlate all position reports 
obtained from the latest sweep of an area with all previously uncorre- 
lated position reports and all previously established tracks in order 
to determine all possible track extensions. A track reliability factor 
is then assigned to each extension to express the probability that this 
extension is the correct one. In this model track extension is primarily 
a function of the selection rule of velocity consistency as discussed in 
Section 2 while the track reliability factor includes those factors 
which must be considered in making a probability determination. What 
these factors are and how they enter the computation of the reliability 
factor form the logic of the probabilistic model. In the discussion to 
follow, the notation developed in Section 3 will also apply, i.e., ia 
is the set of all position reports, Pi ds Di eaOOE Cr Pin reported from 
the it? sweep. 

In order to determine the reliability factor used in this model 
the possible results of correlating any one position report into a track 


extension are examined. For any point P; j from the set P, _ and a 
M 9 


current track A, the set of all possible events in extending track A to 


the set P, is: 
i,n 


Bie Py j is the true extension of track A, or 
9 


[| 


E.. 2 


12 Py j is not the true extension of track A. 
9 


The probability of the event E) occurring, P(E,), is 


P(E,) = 1- P(E, ) 
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Upon further examination the event E, can be seen to include the 


following mutually exclusive events: 


Eo: Pi j is the true extension of some other track X, 
b 
while some other point P. is the true extension 
3 


of track A. 


B.: Ps j is the true extension of some other track X 
3 


and track A has terminated since last sweep. 


LZ? Py j represents a ship which has put to sea since 
b 
the last sweep, with no previous track, and some 


other point Ps y is the true extension of track A. 
b 


By: Ps j has no previous track and track A has terminated. 
9 


The probability of E, occurring is therefore 


P(E 


) P(E.) + P(E,) + P(E) + P(E.) 


and 


P(E) ) Les i P(E.) + P(E3) ~- P(E,) + P(Es) ) (Z,~<11) 


It should be noted that inherent in equation (4-1) is the assumption 
that all ships are detected and their positions reported on each sweep 
by the surveillance satellite. 

The value of P(E, ) obtained in equation (4-1) should be the value 
of the reliability factor used in this model. However, none of the 
probabilities in the right side of equation (4-1) can be computed directly. 
The value of P(E, ) can at best be estimated. In the sub-sections to 
follow, equation (4-1) serves as a guide for obtaining a best estimate of 


P(E, ). 
at 





4.1 Probability Based on Velocity Consistency 

This section sets forth the methodology for determining a probabil- 
ity of correlation based on the selection rule of constant course and 
speed. Basic to this methodology is the assumption that course change 
and ved change are independent events. This assumption is based on 
the observations made in section 2.1 concerning circumstances governing 
course and speed changes. 

If the probability density function for speed change is continuous 
the probability that a ship will change speed exactly k knots is zero, 
regardless of the value of k. In order to obtain a positive probability, 
a small margin of error, /\ x, must be introduced. If tes (sc) is the 
probability density function for speed change, the probability that a 
ship will change speed by an amount k - Ax is 

k+/\x 
f(x) dx 
k=/\x 
However, this computation will still yield a very small probability 
even when the speed change, k, is zero, the most likely value. In 
order to determine a more meaningful value, therefore, the probability 
of a correlation p,, showing a speed change of k +/\x as compared to 
the probability of a correlation which shows the most likely speed 
change, 0 fap is used. 
This probability of correlation is computed as the ratio of 


these probabilities, l.e., 
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k+/\x 


f ,(x)dx 


k-/\x 


*s : +/\x 


fs (x)dx 


- /\x 


However, for constant ae in both numerator and denominator, this ratio 


can be assumed to be; 


= f,(k) (4.2) 


= £80) 


The probability of the correlation based on course change 
Pg , is developed in precisely the same way and if £ A(x) is the probabil- 


ity density function for course change then 


; Io (k) 
Ge” Ta (4.3) 


Ig (0) 
where k is some representative value of course change. 

The probabilities, p, and Po » represent two independent estimates 
of the probability for a valid track between position reports as a func- 
tion of consistent velocity. A best estimate of this probability is 
some linear combination of the two. In this model it is assumed that 


this best estimate for a given course and speed change is P, where 


6’ 


Ps + Po 
> ae 
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The development of Ps Q thus far applies only for those tracks 
that previously existed in the model for which a course and speed 
were known. New tracks may be initiated in the model by connecting 
position reports from the previous sweep of the area which cannot be 
correlated with new position reports from the current sweep. In such 
a caSe there is no track history from which to determine a value of 
course and speed change. The discussion in section 2.1 implies that 
a truncated normal probability density function exists for ships! 
cruising speeds. For a correlation between a new position report and a 
previously uncorrelated report, the probability, p,, is used in this 


model instead of PsQ » where 


£, (k) 
fo (fle) 


and cc) is assumed to be the truncated normal distribution with 


Po 


[ke = 14 knots and Gg. = 2.5 knots and where k is the speed indicated 
by the time and distance between the position reports. This is to say 
that in the track extension phaSe, 1.e., the phase for extending 
established tracks, PsQ is assigned to all possible extensions in an 
effort to determine the most likely, while in the track initiation phase, 
P, is assigned to all possible correlations between new position reports 
and previously uncorrelated positions. 

For the evaluation of this model the density functions, f,(x) and 
fg (x), were assumed to be truncated normal with the parameter values 
Eos = 0, Ge = 2.5 for s in knots and Ho = 0: Og = 2 Lor O in degrees. 
No statistical justification is offered for the use of these distributions. 


They were Selected because they seemed to be reasonable in the experience 
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of the authors, The methodology, however, would be the same re- 
gardless of the distributions used. It may be that experimentation 

or thorough statistical analysis will yield distributions which provide 
more accurate likelihood. These distributions, however, appear to 
adequately represent the conditions of velocity consistency for ocean 
going ships. 


Using these a normal distributions the desired probabilities 


ve: e/b| 
- ae eels 
7 


rqer (02) 


and ina similar fashion 


Pa = eAlg|* 


1 eA 


where 9 = the amount of course change. 


(eds 3) 





(ae) 


k = the amount of speed change. 


s 
" 


the speed indicated between position reports. 


and 


t 


the truncation factor, 
4elel Probability Based on Mutually Exclusive Events 

The factor P59 expresses a likelihood that a particular track 
extension is a valid one based on the consistency of course and speed. 
In the probability sense, the set of all Ps g for the extension of some 
track A cannot be considered as a mass function to determine which 


extension of A has the highest probability of being valid. This is the 


EB 





case because, in general, the set of all p,g for track A will not 
sum to unity. However this situation may be overcome by assuming that 


all possible extensions of track A from Pima to P and the associated 
3 


alee 
Psg form a set of mutually exclusive, exhaustive events. Thus a likeli- 
hood factor that represents the probability of a track extension being 


valid in view of the k possible extensions can be computed as 


LAP: 42, 3) a (PsQ ) 
Pa FS 9 Fe 
where 

LAPS 43 5) is the probability of Pay j being the valid 
extension of track A with respect to consistent 
velocity and all other possible extensions of 
track A. 

P59 ) is computed as in section 4.1. 

k is the total number of possible extensions of 
track A. 


It should be noted that the relative magnitude of the Ps have 

been preserved in the construction of L(Ap, 44 j For example: if 
9 

k= 1, L(Ap, ) = 1.0 and if k = 3 and these three possible track 


ats all al 
extensions are equally likely then L(Ap, ) = .333 forall thiee 


se ee 
regardless of the value of Pso ° 
This same development applies to the computation of the probability 
of a newly initiated track being a valid one based on the speed indi- 
cated and weighted by the number of other possible new position reports 
which could be correlated with the previously uncorrelated position. 


The only difference is the use of p, instead of Pso wherever the latter 


appears 6 
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401.2 Probability Based on Cross Mutual Exclusion 
The factor L(Ap, ,4 3) is developed by a consideration of all 
3 


position reports Ps . from the set Poa 


1. j which fall within the dis- 
bf 


1, 


tance Dake of the last known position of track A. It is also possible 


that the position report P; is considered as a possible extension of 


ree) 


more than one track, say tracks A, B, and C. If this is the case then 
. there is an associated Ps as computed in section 4.1 for each of these 
possible tracks. The arguments presented in Section 4.2 for the develop- 


ment of L(Ap. 


.) will again hold and a new probability factor L(p. 
ate ee i 


ee 


can be defined which reflects the velocity consistency of Ps with 


olla 


respect to track A normalized over the other tracks for which P; could 


te 


be the extension, i.e., 


wnrFrn 


1) ey, 
eo ae Ty, J 


hs 


(Pg ). 


all 


where 
L(pPi 41, 54) is the probability of Pidl, j being the valid 
extension of track A based on consistent velocity 
and all other possible tracks for which p, 
oR SUG) 


could be the valid extension. 


PQ is computed as in section 4.1, and 
ig is the total number of track extensions containing 
Pi+l,J° 


t thi int the difference between LiAp. and L(p: A 
At this poin ifference betwe ( Pign, 5? (Psd, j ) 
should be noted. L(Ap, 15 ;) is the probability that position report 
is the valid extension of track A with respect to all other 


Seo) 
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position reports that may be valid extensions of track A. LPs 44 j4) 
b 
however is the probability that position report Piy j is the valid 
3 
extension of track A with respect to all other tracks that may contain 


Pi+l,j° 


The following computations can now be made: 


L(Ap: 


i41, 5 the probability that some position report, 


Say Pit] ke other than position report Pi+l,j 
is the valid extension of track A, 
= 1- L(Api 41, 5); 


and 


L(p, 44 jA) the probability that position report pj 4) ; 
9 9 
is the real extension of some track other 


than track A 


1 - L(pi 41, j4)- 


Referring to equation (4-1), P(E,) can be computed as: 
P(E,) = MAPS aj? x EN st pe 


42 The Probability of Track Termination and Initiation 
Neglecting the last three terms of equation (4-1), a good estimate 
of P(E, ) based on velocity consistency and mutual exclusion with respect 


to both track A and Pray j is given by 
? 


P(E 1 - P(E,) 


_ = 


This approximation of P(E, ) does not include the possibility that track 


A has terminated or the possibility that p,__ j represents a ship which 
3 
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has just put to sea. Therefore, the possibility that tracks may be 
initiated or terminated between sweeps must be admitted for every 
track at each track extension phase, 

Probability of track initiation. If A is the mean rate of 
leaving port for ships all over the world then 1/) is the mean inter- 
val between ships leaving port anywhere in the world. The average 
number of ships leaving port in a time interval t is therefore ty ° 
If N is the average number of ships at sea the probability that any 


particular one of them left port in the time interval t is given by 


t/) t 


MnPiia,j? =a NX 








where, if t is the time interval between sweeps, then GS sie is 
the probability that the ship represented by the position report 
Pi4l, j put to sea between sweeps i and i+l. The justification for 
computing LPs 473) in the manner shown above is the observation 
that ships put to sea in completely random fashion. This assumption 
is false at any particular port because the rate of departure depends 
on conditions of weather, tide, and daylight or darkness. However, 
over the world as a whole the effect of these factors at one port is 
assumed to be cancelled by opposite conditions at another. There are 
not an infinite number of ships as would be implied by a completely 
random distribution; however, there is a relatively large number of 
ships in existence compared to the number of ships at sea, and the 
turn around time as compared to voyage duration is relatively short, 
The effect i similar to an infinite supply of customers in a queue. 


“In the evaluation of this model the parameter values of A and N 
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were chosen as .10 and 5000 respectively. 
Probability of track termination. In order to determine the 
probability of a track terminating between sweeps the authors assumed 
a truncated normal distribution to represent the distribution of 
voyage lengths of ships at sea. Using this distribution the probability 


that a track will terminate between sweeps i and i+l, L(-A) can be conm- 


puted as 
d 5 
U-a) = Let ae 
oS SS t | dx 





where d is the total distance for track A from the point of initia- 
tion up to its position on sweep i+l. 

There may be some justification in claiming that the probability 
of termination should be conditioned on the fact that the track was 
extended on a previous Sweep, i.e., that the probability should be the 
probability that track A has terminated before sweep i+l given that the 
track A had not terminated before sweep i. This conditional probabil- 
ity was not used in this model because the premise that the track had 
not terminated before sweep i is not a known fact. In the framework of 
the probabilistic model the track's existence on the last sweep is only 
a probability. 

Based purely on the experience of the authors, the parameter values 
of mean voyage length jin = 3000 nautical miles and 0, = 1000 nautical 
miles were chosen to test this model. More accurate parameters could 
be derived from a careful analysis of shipping records. 

It should be noted that if this model were to be expanded to include 


a computation of the distance of a position from known sea ports then 
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both L(-A) and Lip, 4) ;) would be zero if the distance from port 
9 
were found to be greater than Sax? 
The values of P(E3), P(E, ), and P(E.) in equation (4-1) can now 


be found by using the values of L(A) and Létp; 44 js i.en, 
) 9 


— 


32> mst 


P(E) 


o> 
As 


tt 


L(pi 41, 34) x L(-A) 
P(E,) = L(Aps 4) 5) * L(-pyyy 4) 


P(E.) = L=Pi yy 5) x L(-A) 


4.3 The Track Reliability 
Utilizing the individual probabilities calculated in the preceding 
sub-Sections, the complete formula used in estimating the probability 


that a particular correlation is a ‘valid one, P(E, ), is 


P(E, ) =1<= L(A 43 5) x LPs 47 54) o Mpa ya yA) x L(A) 


This estimate is the track reliability factor used in this model to 
determine both the most likely new tracks and the most likely extensions 
for established tracks. In the next section a detailed sample problem 
will be used to demonstrate how track reliability is employed. 
4e4 Application of the Model 

The probabilistic model was programed in FORTRAN 63 to be run on 
the CDC 1604 computer. Summary flow charts and a program listing are 
included in Appendices III and IV. The logic of the program can best 


be shown using the set of diagrams in figures 6 through 11. 





For illustration purposes assume that in the initial orbit of 
the satellite, the set of position reports aL is obtained. Since 
there are no previous points with which this set can be correlated, 
they are retained in the model as uncorrelated positions. After the 
second sweep of the area a new set of positions 2,3 is recorded, It 
is determined that each of the points in the set Po 3 is within the 


radius S,, of each of the previously uncorrelated points in set Pi 3° 


Figure 6 shows the nine possible tracks which must be considered. 
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Initial Tentative Tracks 


Figure 6 


A value of P(E) is computed for each of the nine possible tracks 
and all nine are retained until the next sweep. 

At the third sweep a new set of position reports 53 are received 
and the selection rule of consistent velocity can be applied. As before, 
all positions from the set ays which fall within the Snx radius of the 


positions in P, 3 are considered possible extensions of the nine tracks 
9 
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established in sweep 2, and P(E) is computed for each. In this track 


extension phase the probability Psg can be used instead of P, which 


was used in the previous phase. 


3,2 


3,3 





Most Likely Tracks 


Figure 7 


In Figure 6 the most probable tracks are indicated by the solid 
lines. Utilizing the selection rule of consistent velocity to extend 
tracks the correlation (p p p, ,) has a higher track reliability 

Bim ol mS) 3) 2 
factor than any other combination involving p, ,. Hence (p p ) 
1,3 1,3 P2,1 
and (P) 3 P 9) are discarded. By the same reasoning (Py 4 Poy) 
(Py 4 Po 3) (Py 9 Po os and (Py 2 Pa 3) are eliminated, Only three 
of the tracks produced by the position reports from sweeps 1 and 2 
remain. 
The possibility that (p, 1 Po >) represents an initial track with 
9 9 
a course and speed change to produce the track (py p, ,P, _) must also 
A” 2c cn oes 


be considered when extending tracks from PL 3 to P; 3° The new possible 
9 


9 
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track extensicns from P t u 0 

s from lies O tee to sii due to "near" constant 
velocity alone are illustrated in figure 8 by solid lines while those 
tracks possible due to a speed and course change are illustrated by a 


combination of solid and dashed lines, 





Tentative Extensions of Most Likely Tracks 


Figure 8 


For the tracks in figure 8 the solid lines receive a higher 
reliability factor than do the solid = dashed combinations. Figure 9 
illustrates the next phase of extension where the possible tracks are 


correlated with the elements from From the new set P, 3 the 
9 


Phat 


probabilistic model has assigned only two tentative track extensions 


which will be used in the next phase of track extension, 
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4,1 
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Firm Track Established 


Figure 9 


In the foregoing example tracks were initiated by the correlation 
of position reports from the current set Pi+lin with uncorrelated 
positions from the set Pi me This type of track initiation will occur 
in the first few orbits of the surveillance satellite and may possibly 
occur at any time later. A much more likely situation in later orbits 
is that in which the position from the new set Pitlin which represents 


a ship which has just put to sea will be considered as a possible ex- 


tension for one or more established tracks, It would not then be carried 


over to the next sweep as an uncorrelated position report. 
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Track Split 


Figure 10 


In figure 10 two tracks are shown,each with three possible ex- 
tensions, One of the points from the set Pielin must be a position 
report for a ship which has just put to sea. The decision as to 
which of the positions represents the new track is postponed until 
the next sweep. Figure 11 includes the position reports from the set 


P After computation of P(E, ) for each of the possible extensions 


ego : 


track 1.1 extended to Piao] is determined to be the most likely 
2 
extension of track 1, and track 2.3 extended to p, |, 3 the most likely 
9 


extension of 2. This leaves the correlation (p ) as the 


i+1,2, °1+2,2 


new track, 
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New Track Initiated by Track Splitting 


Figure 11 


The track correlation problem for the probabilistic model as 
illustrated in figures 6 through 11 is, of course, for the particular 
sets, Pion’ However consideration of these figures should give a 
basic understanding of the model. 

The correlation process is a matter of examining three consecutive 
sets of position reports to select a most likely combination of posi- 
tions. All but the most likely combinations involving the earliest 
position reports are eliminated. However, in going from one track ex- 
tension phase to another all combinations involving the last two sets 
of position reports are retained in the model. It should be noted that 
the elimination of the least likely backward extension is the equivalent 
of establishing the track's identity with respect to tracks previously 
established. If a track cannot be identified with a previous track, it 
must be assigned a new designation and considered to be a ship which has 
just come under surveillance. It should also be noted that retention of 


only the most recent position report together with a velocity vector is 
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the equivalent of retaining the last two position reports and a pre- 
dicted position for the next sweep. Thus, an identification symbol, 
a last known position, and a velocity vector provide the necessary 
and sufficient information for tracking in this model. 

4.25 Tests and Results 

In order to test this model a CDC FORTRAN 63 program labeled 
TRACK IV was constructed to conform to the logic of the model. This 
program and associated flow charts are included in Appendices III and 
IV. Another FORTRAN 63 program was developed in order to supply initial 
test data to TRACK IV. This input program is called READOUT and is in- 
cluded in Appendix III. 

The only formal tests imposed on the model involved variations on 
the situation shown in figure 12. While the results of these tests 
provide a basis for qualitative observations there was insufficient 
time to impose the variety and quantity of problems required for a 
quantitative analysis. 

The situation illustrated by figure 12 involved three tracks con- 
verging to an area and changing course when in close proximity. At 
about the time of course change for the three original tracks two new 
tracks were introduced. This problem was posed to TRACK IV on ten 
different runs with the interval between sweeps varying from one to 
twelve hours. Tabulated results of these runs are shown in table 13. 

"Time to Resolve" is measured from time 64 when the problem became 
critical until the time when all tracks had been firmly identified and 
all course changes had been properly recorded. 

The column entitled "Tracks Correctly Extended All Orbits" applies 


to tracks properly initiated and correctly tracked without loss of 
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identity throughout the whole of the problem. Improper tentative 
extensions may appear on any particular sweep with a higher track re- 
liability than the true extension. However, if on the next sweep the 
true extension is assigned the highest track reliability the track is 
considered to be correctly extended. It should be noted that the most 
common error which precluded a track from being listed as correctly 
extended was that at some point the track was redesignated. In such 
cases the track was never really lost for more than one orbit. Cases 
of redesignated tracks are tabulated under the next column, "Changed 
Identity", 

An entry is made under "Changed Identity" each time such error 
occurred. In the entries for the two hour interval one track was re- 
designated three times and another track was redesignated once. All 
of the errors which occur in the three hour interval concerned only 
one track. 

The column entitled "Improper Extensions" refers to cases in which 
a proper track jumped to assume the identity of some other valid track 
while the original track was either lost or assumed a new identity. 
This is considered to be the most serious type of error. 

The column heading "Time to Resolve New Tracks" refers to the two 
tracks introduced midway in the problem, The time is measured from the 
first appearance of the positions representing them until TRACK IV holds 
them as identified tracks with the proper extension being assigned the 
highest reliability. 

The last column lists those occasions where track designations 
were assigned to a group of possible correlations none of which were 


correct. On all three of the occurrences listed, the false tracks were 





dropped on the next sweep. 

Another type of error was discovered in the tests which is not 
tabulated. In the one and two hour sweep interval runs,a number of 
tracks appeared which were duplicates of other firm tracks. This is 
due to a shortcoming of the program, TRACK IV, and is not a consequence 
of the logic of the model. The shortcoming is easily rectified by the 
addition of a sequence in the program which searches for duplicate 
tracks and eliminates the least likely. 
4e6 Recommendations and Conclusions 

The tests conducted indicate that the probabilistic model can be 
developed into a very useful instrument for assisting in the tracking 
of ships on the high seas. Only two of the shortcomings tabulated 
are inherent in the program logic. These are improper track extensions 
and the originating of false tracks. The latter is not considered 
serious because in each case the false track was in the system for only 
one sweep. Improper track extension occurs when at a particular read- 
out the inter-relationship between position reports indicates that two 
ships are making improbable course and speed changes. This fault can 
never be completely eliminated from the model but implementing some of 
the recommendations to follow may greatly reduce the frequency of 
occurrence. 

As a result of the tests conducted the following modifications 
should be made to TRACK IV before further tests are undertaken. 

1) More sophisticated programing methods should be employed in 
constructing the correlation lists. In the present model a maximum 
of 29 correlations can be tabulated for each track. The number of 


possible correlations increases as the square of the number of tracks 


48 





in an area of congestion. Thus in its present form TRACK IV reaches 
a saturation level when 5 tracks CONVES Ze o «. 

2) An internal track strength indicator should be assigned to 
each track. This can be a simple tally of the number of times a track 
ls successfully extended. The track strength indicator will be used 
in connection with the next two recommendations. 

3) In its present form TRACK IV computes course change and speed 
change by taking the difference between those indicated by the last two 
positions and the last course and speed. This is accurate only if the 
ship changed course at precisely the time of the last position report, 
a very unlikely event. The result is that any course change appears as 
a Simultaneous course and speed change neither of which are correct. 
This reduces the track reliability to a level far below what it should 
be, and probably is the cause of most cases of changed identity and 
improper extensions. A routine should be added to TRACK IV so that 
when track strength indicates that the track has been successfully ex- 
tended a minimum of perhaps three times and a course change is indicated 
with a simultaneous reduction in speed, the speed wiil be assumed to 
have remained constant and the exact amoint of the course change computed 
by solving a triangle with the law of cosines. 

4) Identity changes occur because of a shortcoming in that part 
of TRACK IV which eliminates the least likely backward extension of 
tracks previously listed as tentative. It is conceivable that all 
possible extensions for a track appear with a higher reliability as 
possible extensions of other tracks. In this case all the possible 
track extensions carrying the proper identity are eliminated, and the 


proper extension of that track appears as a split for some other track. 





It is then handled as a new track. This situation can be remedied by 
inserting a change in the sequence which handles track splitting so 
that the newly initiated track is compared with all tracks which have 
been dropped to see if it is in fact a continuation of one of them. 
If so and if the strength indicator of the old track was greater than 
three then the old track's identity and history should be adopted for 
the new one. 

5) Track duplication occurs as a result of several procedures 
in TRACK IV. The best way to rectify the situation is to add a sequence 
at the point in the program when all new track histories have been 
prepared. This routine would check for duplicate firm tracks and retain 
only those with the highest track reliability or those with the highest 
strength indicator. 

With the implementation of these recommendations the probabilistic 
model will be capable of handling much larger and more complicated 
problems and will yield much more reliable performance. 

The most surprising result from the Series of tests of TRACK IV 
is the effect of extending the sweep interval. If those items used as 
headings in table 13 are used as criteria for performance then the ten 
hour sweep interval appears to be the most efficient from most stand- 
points. However, the results of the runs for sweep intervals between 
two hours and eight hours indicate that the factor which affects 
performance is the particular relationship existing among position 
reports on successive Sweeps rather than length of the sweep interval 
alone. The model will have to be tested with a large number of varying 
problems before any significant observations with respect to sweep 


interval can be made. 
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5. Recommendations and Conclusions 

The following recommenda tions and conclusions are based upon the 
analysis of the evaluation of the proposed models. 
5.1 Recommendations 

The validity of the correlation models appears to warrant their 
further expansion and investigation to determine: 

1) The effect of relaxation of the assumption that all ships on 
the ocean and within range of the sensor will be detected and reported. 
In particular, further study should attempt to ascertain the lowest 
blip-scan ratio that will permit an acceptable level of correct corre- 
lations. 

2) The optimum size of the areas constructed about the track 
projections. Statistical analysis based on many random simulations 
should permit this determination. 

3) The maximum shipping density (number of ships per unit area) 
that will permit an acceptable level of correct correlations under 
various constraints regarding the number of course and speed changes and 
orbit interval. 

4) The effect upon the maximum shipping density of the input of 
externally gathered information concerning the course and speed of a 
specific ship. 

5) Whether these models can accept data gathered and preprocessed 
in the manner of the Pacific Missile Range simulation and correlate the 
data correctly. 

6) The appropriate measures of effectiveness to be applied to this 
model to facilitate analysis of the meaning and effect of various para- 


meter changes. 


at 





7) The design and effectiveness of a new model combining the 
essential features of the models proposed in this thesis. 
5.2 Conclusion 

The results of the evaluation of the two models proposed in this 
thesis indicate that track correlation using satellite gathered infor- 
mation can be satisfactorily accomplished by these and by similar 
correlation models, It is anticipated that further refinement of these 
models would provide a track correlation logic with very high capability 
and would be of significant value in further studies of the feasibility 


of ocean surveillance by satellites. 
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APPENDIX I 
FORTRAN 63 PROGRAM-DETERMINISTIC MODEL 

This appendix contains the computer program which was written to 
permit computer simulation and evaluation of the deterministic model. 
The program is written in the Control Data Corporation FORTRAN 63 
computer language and the model was evaluated by computer simulations 
on the CDC 1604 computer. An example of the problems generated by the 
program and a discussion of the input and output associated with the 
program are contained in Appendix II. 

Comment cards are used extensively in the program to explain the 
flow of the problem and to define the variables used. Detailed flow 
charts of the program logic are on file and may be obtained from the 
Chairman, Department of Operations Analysis, U.S. Naval Postgraduate 
School. 


The program uses the following FORTRAN 63 library functions: 


1) RANF 

2) SINF 

3) COSF : 
4) ATANF 

5) SQRTF 


The program compilation time is 2 minutes and 53 seconds. Com- 


putation time ranged from 8 - 18 seconds per simulation. 





*SNOTLISOT GSLVISYYOOINN ONINIVINOD @ 30 SMOY JO YAEWNN=YSSYNNN 
*GSLVNOISSG SAOVYLl JO YAEGWNN=ADVYLN 
°F SWIL LV I Adval SV AGSLVNOISSGA dIHS 3O YSEWNN=(FSITIADVYL 
°") SONVHD G33dS JO LNNOWV=(1)5S9VW 
°) SONVHD 3SYNOD JO LNNOWV =(17)D59VW 
°J0VW 38 OL SSONVHD G33dS JO YSEWNN =QdSON 
°"30VW 38 OL SSONVHD 3SYNOD 3O YSEWNN =SNDON 
*SJONVHD G33dS HL-X SHL ONIDVW dIHS 3O YSEWNN =(2X)SDSWVN 
°SJONVHD 3SYNOD HL-A SHL ONIAVW dIHS JO YSEWNN =() DD3SWVN 
*YOLVYANIAO YAGWNN WOGNVY TIWWYON V SI E9AZQNY 
*YOLVYSNASD YAGWNN WOOGNVY WYOSINN Vv SI S4NVY 
°SNILNOY ONTLVYSNS9 YSEGWNN WOONVY NI G3SN 3YV SINAN ONV XW 
*G3SLVYSNS9 38 OL SW3180¥d WOGNVY JO YSEWNN=GOUdN 
*dIHS ANV Q3L1LIWYSd G35ddS WNWINIW=QdSNIW 
y *dIHS ANV G3LLIWYSd G3S3ddS WNWIXVW=ddSXVW 
*SSJONVHD Ga3adS 3O NOILVIASO AYVONVLS=A3500dS 
*SSONVHD 3SYNOD JO NOTLVIASZG GYVUNVLS=A3Z0SND 
GIHS 17V 3O NOTLISOd IWWILINI 3SHL ONINIVINOD VSuvV SYVNOS SHL JO 3ZI1S=0S3ZI1S 
°GQ3SVSYONI N33EG SVH 3ZIS V3eV SSWIL 30 YSEWNN SSLVDIGNI S3TDADN 
8 - JO + SV HLIGOIM V39V INIOSTYO SSHSIIGVISS G-IdOML=ODHLGIME GNV G=OHLOIMG 
°“SLONA V - ¥O + SV 3ZIS V3AYV AWNIODTYO SHL SSHSINEGVLSS V=OH1Ld508 
°03SS390¥d 38 OL SSSVHd JO YSEWNN=VON 
°03SS3908d 38 OL SLIGYO JO YAEGWNN=LIGYON 
*SWIL ONILYOd3SY HDV3 LV G3LYOdSY SdIHS JO YASWNN=LISOdN 
*SSONVHD G35dS ANY 3SYNOD N3SML13E SIVAYSLNI SWI 
JN3S¥3SSTQ HLIM LIT 1931 OL WVYDOYd AHL ONIGN3SLX3S 4O NOIIVdIDILNV 
NI QSQNIDNI SI SIHL °*S1LYOd3yu¥ NOILISOd N33SML35EG JWIL=(1)LV1L730 
*SLYOd3Y¥ NOTLISOd NS3SM1598 3WIL=3wIl 
DISOVWSDDSOVWSAdGSNIWSAOdSXVW W354 SJdAl 
1*HDVYLl YSOSZLINI SdAl 
DOINANS LISOdGNS GS DHLGIMET 
*HIGIMESH1Ld308* YSSYNNN SOV ISNS VANS SLIUTSONS XDVYLSADVYLN NOWWOD 
(SSO0TISNDS COTIDSOVWS (OT) DIOVWS (OT) SDISWVYNS (OT )DDIWYNZ 
pee lee OT) Oe OL) IVIingasis*S2iyovuls (9C) 14 (e*S*Olivi 
“(9TSOTSOTIGS(SSOTIAS S(SSOTIXS £(6*40T)0dSS*(SSOTISNIS NOTSNIWIA 
, ONTIAIVAL WVY9DOdd 


_— —_ - -_— ~ - - ~ - - — - —_ ~_~_ — - - _ * - _, ~ - - - - t - 


VUUUVYUUUYUYUUUVYGUUUYUUUUUVUUVUUVUUUUUY 





°4+ OQNV £2 SASVHd YOS SAITOH LVWYOS YVIIWIS 
*SNOILISOd Q3LVI3SYYOD YOS NOTLVYNOIS3SO ADVYL SNIVINOD (#*rsI)8 
*T+N LIGYO WOYS CF NOITLISOd OL N LIGYO 
WOYS I NOILISOd WOUdS 1359 OL AAYINOAY ASYNOD AHL SNIVINOD (Eesrfs1)a 
Shwe hao woods ef  NOwemISOd Ol N Iledo 
WOY¥S I NOILISOd WOYNS L359 OL AGSYINOAY AAAdS SHL SNIVINOD (2fr*I)e 
°YSEWNN SASVHd SHL SNIVINOD (TSFS1)8 
Daher erOewOds shdOdsad NOTEISOd SHLl 30 YaagWNN SHI OL SdYssayu £ 
*N LI8yO WOYUS LYOd3SyY¥ NOLISOd S3SH1 3JO YSEWNN SH1L OL SYASAY 1 
°% 3SVHd YOS NOILVWYOSANI SNIVINOD OTSET=X 
°€ ASVHd YOS NOILVWYOSNI SNIVLNOD Z2T*6=% 
°c S3SVHd YOS NOILVWYOSNI SNIVINOD 8*S=%H 
*T S3SVHd YOS NOITIVWYOANI SNIVINOD 7fT=H 
--NOITLVWYOSNI ONIMO11IOS SHL SNIVINOD (XSI) 
eS) SIA ee hel 8, 
WOdamasOnml Ol 2 L19@¥0 WOUS SLYOdaY NOILISOd N3SAM1S8 SNOTLV IAN SNYSINOD 
€ ASVHd °*2 LIGYO WOYS ASOHL OL T LIGYO WOYS SLYOd3aY NOILISOd N3AM1L3¢E 
SNOTLV17S¥Y¥ SHL SNYSDNOD T ASVHd SYSHM SASVWHd AG ASSSADOYUd SI NOITLVWYOSNI 
"dNlasS SIHL YSONN AONVYDISINOIS ON SAWH (TS1)0dSS INV (T*I)SNDS--3LON 
ef -EL edo we 
JWIl 1V I dIHS 3O SLVYNIQHYOOD A ONVW X BSHL SLVDIGNI (FSI)AS ONV (FST) XS 
°f OQNV T-f LIGYO N3Z3SM138 
AWIL ONIYNG I dIHS 3O Q3a3dS ANV S3SYNOD SLVDIGNI (F*S1IGdSS GNV (fF*1I)SNDS 


°"O1S30Q ANILNOYBNS AG O 

JWOS WOUS a3Aowae Noo deo NAN Is aNOeISVan IV IVHI SSLVYSIGNI [=ala7a0N 
°g@ JO NWN10D YO MOY N3AIO V NI SLNSWS73 JO YSEWON SHL=W373N 

*LN3W3A173 3NO A INO 

SVH HDIHM O L3S WY WOYS SI GSLDSCOYNd ONISS ADvYLl SHL SASLVDIGNI T=S0INNAN 
*“QSNIWVX3 ONISG 

VaudV SHL NI SI LNAWS3173S G3LVNOIS30 SNO LSV37 LV SSLVDIONI T °39* SVTSANN 
*VAYV NSAIO WV NI dIHS ASANO NVHL SYOW SALVDIOGNI T=9V14SF 

*MOY YO NWN1IOD DJISIDAdS V JO ADAHD-NI GSIdNDDO SV3SYV ON SSI TdWI O=NI971 
°V5S9dV NSAID V NI SdIHS 3O YABWNN=NIGN 


YVUUUUVYUUUVUUUUVUUUVUYUUUYUUYUUUYUUUYUYUUUUUUUUUUUY 


oA 





; SNNILNOD 
IDOMLe¥SNIV=E(VFSVIISNDS $F$ LIQGYONST=VF 8OT OA 
(T-)ANVY=SNIV $F LISOdNST=VI SOT OG 


*IdOML OL O WOYS AIWYOSINA 3SYNOD WILINI 31VY¥3N39 


SANTINOD- 


O° O002+tOSSZISe(T-) ANVYU=H(TSVIIDASSO°OO0Z+OS3SZISx(T-) ANVU=(TSYVI) XS 
LISOdDN*T=VI LOT Od 
OT+JSINQAN=S3INNAN 

Fi XWST=I TST OQ 

(T-)ISNVY=S 

XW*ST=I OST OO 

G¢+XW= XW 


ec 7 ls =3WILH9*xz*z7°9345 =OSAZISHBSXZ*2°SGAS =AAIGGdSHBS XZ2SET 


°94* =AZ0SNDHBSXC*7°9S3% =AdSNIWHBSXZ*7°94* =AdSXVWHE* XOT*//) LVWYOS 
IWILSOS3AZISSAZ0GdS S$ AAGSNDSAdSNIWSAdSXVWS2246 INIUd 

(E€1S SOHESXES ETS YISWAN W31808d WOONVYHZZ*XOTS///* THT) LYWYOS 
BOYdNSNNYWANS VET INIYd $ POYdNST=NNYWNN 2ST Od 


°OS3AZIS *31GVIYVA SI 3YVNDS 3O 3ZIS *O0Z=A * OOZ=X Lv 
NISTYO HLIM JYVNOS NI GalNGIYLSIG AIWYOSINA NOTLISOd IILINI 31VYe3N39 


*SNOILISOd IVILINI GONV *S033dS *S3SYNOD Y3Z1N3 


T-VAN=ON $ T-LIGYON=VAN 


€6T*°9=ODHLOIMAS60* O=OHLOIMESO® 2=OH1Ld309 © 


O° S=QdSNIWSO°SE=AdSXVWS0°S=AZ00dSS$S °O0=AZIOSNISO*OOT=0S3Z1S 
OC=GOUdNST=SINOANSLTI=XWSO°T=3WI LSSG=LIGUYONSOT=LISOdN 


*SYaLaWVuVd YBiN3 
GIS6ESB*L=VId $ ESBTEBZ*9=IdOML $ SZ6SGTHT*E=Id 


SINVISNOD Y31N3 


BOT 


Zot 


OOOT 
TST 


OST 
Tool 
C276 


vEeT 
esl 


UUYU 


VUUUYUY 


UYU UUYU 


57 





JNNILNOD 

| IJNNILNOD 

OZT*ZZI ( (4) SDIWYN°OS* (VX)SDSWVN) JI 
TZt*OZT (x*oO3*vy) SI 

4*T=Vy Oz2T oa 

O° T+0°OTx(T—) INVY=(X)SDIWVN 
QdSDN*I=xX 611 Oa 

IJNNILNOD 

3JNNILNOD 

STISZIT( (9) DDSWVYN°O 3° (VX) DDSIWYN) SI 
9II*SII (x°03*° VX) SI 

X*L=vy¥ SIT oa 

O° T+O*OT*(T—) SINVY=(A) DDIWVN 
SNDDNST=xA +I OG 

O° Z+0%L*(I-) INVY=GdSIN 

O° Z+0%L*(1T-) 4INVY=SNION 
3NNILNOD 

O-=(4)SD3WVN 

O-=(4)D39WYN 

LISOdN‘T=X 8TI OG 
LIGUON*‘E=F EIT OG 


*SOSY3HL SQNLINOVW 3SHL GNV SSDNVHD 3H JXVW 


SdIHS HDIHM JNIWYNSL390 N3HL °((8*Z)WYOSINN) 3QvW 39 OL SSONVHD G33adS GNV 
3SYNOD 4O Y3SGWAN 3BHL JLVYSNSOD Z LIGUO OL LNINDISENS GOIY3d AJWIL HDV3 NOS 


“ 


SJNNILNOD 

AZG=(VFSVIIAdSS 
LIGYONST=VF 60T OG 
CTTSOTT(GdSXVW°39°AZ0) SI 
TTT*OTT (GQdSNIW°117°AZ0) SI 
0O°ST+A3S0QdS#ASG=AI3RG 

(ASQS SINAN) EOSAZONY 115d 
LISOdN*T=VI 60T OG 


pest ONY SG 39 GdSs 1VHi HONS *(AZS0GdS*STIWYON WOYUS G335dS IVWILINI 31VYSN39 


6Tl 
Crag | 
Tel 
ccl 
7TI 


STl 
9TT 


ole 


ST 


60T 


ett 
be 


OTT 


VUUUNUYU 


cS 





NHS CS XS SET € 


JANNILNOD 


(€1*S =SSJONVHD G335ddS 340 YIEWNT 


=SIJONVHD 3SHYNOD JO YAGWNNHIO?CSXGSETS =CHE*SXOTS//) 1LVWYOS 


} 


GdSON*SSNDINSCSEET LNIYd 

JINNILNOD 

1° 0-GdSXVWwe (OF *NI)QdSS 

S8ZT*ZET (ddSxXvwe3s9° (OF SNI)GdSS) JI 
8ZIT OL O9 

GdSNIW=(0F*NI)IQGdSS 

TETSOET (AdSNIW*°L1° (OF SNI)GdSS) SI 
(Ol) DSOVWH+(T-F SNI)OdSS=(OF*NI)IGdSS 
LIGYON*‘F=OrF 8ZT OA 

(O11) SD3IWVNENI 

QdSON*T=O!I 82T Oa 

JNNI LNOD 

(OL) 2309VW+(T-P SNIISNDS= (OF SNIISNDS 
LIGYON*r=Or 92T OA 

(O1) DDSWVNENI 

SNDINST=0OI 92T OC 

JNNILNOD 

A3dGdS*AI0=( 1) D9S9VW 
(AJG*SINANIESAZONY WIVD 

QdSON*T=1 S2T Oa 

JINNILNOD 

EZTS42T (1d°19° (1) D59WW) JI 
AZOSNDI*AZOG=( 1) D59VW 

(ASO 4INAN)IESAZSONY WIVD 

SNDONST=1 €2T OG 


ell 


ctl 


8c 
cel 
Tel 


Oel 


Tor 


Sc 


ee ' 


721 


a9 








eset teeta Ti xXOl tt °94* KOS Cee dt xe*cI* LHI*xtl*/)LyWweos To00d 
(LISOdDNS THIS T 
Pradun le Otr lyAS tp * 1 ixs*r* (F< )gdsstir*s1)Snd*1)))*T002 LNIYd 
. (QYOOD AHL*SXOTSal 
YOOD XHLSXOTSAWI LHS XOTSAdSHE £ XOT SSNDHE * XOTS HDVYLHSSXOTS///)ILVWHOS 0002 
000¢ INIUd 


“ S1ISOd “‘SaaadS *S3SYNOD LNO LNIYd 


9oOc°L Gel *1ISHSS=(fF SIISHD O10T 
LIGYONfT=F OTOT OQ 
LISOdN*T=I OTOT OQ 


“ONILNIYd 3YOSSG S$33YN9DTG OL SNVIGVY WOUS 3SYNOD LYZANOD 


(DSIANISH(C ST IGdSSe (IT) LVLTIG+(T-F ST ASHIFSIDAS EOOT 
(DS)ASOD# (CST IGdSSe( 1) LVLITVFO+(T-F 61) XS=(F41)XS SOT 
IdOML-DS=DS ESOT 

WGOTS ESOT CIdOML°399°DS)4I 
(FST ISNIS—-VId=59S TSOT 

€SOT OL O9 
PadOMit= Cl 1SMDS=tF* 1 eshos OSol 
TSOTSOSGOT (IdOML*° 39° FFSEIISNDS)AI 2SOT 

TIC TOL 205 
IdDOML+(FfSTISNIS=(FSIISNDS O9OT 
CSOTSOSOT(OP OP LIP(FSIISADS)ISAI TS9OT 

JWIL=(T) LVL130 

LIGYONS2=r COOT OA 
LISOdN*T=I1 €OOT Od THT. 


SLIGYO VWIv YO4 SdIHS WIV YO4 SLISOd 3LNdWOD 


———-~ 





JNNILNOD 

O-=(FS1T)ADVUL 

Sc*T=I 7c0£ OC 

- G*‘T=f 720E OA 
: O° T-=(vrsvIIQ 
LISOdN*T=VF SS6 OQ 
LISOdN*ST=VI $S6 OG 

O° T-=(VASVFSVIDD 

€*T=VA 4766 Od 
LISOdN*T=VF 7S6 OA 
LISOdN*ST=VI 4766 OG 


°d QNV D 42104832 


0°O0-=(@n*aer*al)a 
ANS T=@X SOOT OQ 
LISOdN*T=8F SOOT Od 
LISOdN*T=8I SOOT OQ 
7x (T-LIGYON) =OAN 


(ASF *1)@ 3210832 


(VF SVIDASH(ESVFSVIIV 
(VF SVIXSH(CSVFSVIDV 

JWI Le VF=aCTSVFSVIIV 
LISOdN*‘T=VI 7OOT OQ 
LIQGYONS T=VF 7O0T OQ 


mise) y ONT StisOd T1iV¥ YSslN3 


72O0L 


SS6 


7S6 


soot 


Loot 


700T 


UUYU 


VUYU 


61 








(THT) LVWYOS EZOE 
((CX2SH7°BSS)OTSXOTS/ILVWHYOS CZOE 
(LISOdDN* T=FS (NSF SI1)8) £220€ INIY¥d TZO0E 

7—- l+Wex7=> 

7*T=1 [TZ0€ OQ 

LISOdN*‘T=I TZO0E€ OA 

ecOe INIYd 

VANS T=W TZ0€ OG 


@ INIUd 


; SJANILNOD 9OOT 
VISHL=(T-FxVS TSX) ASd33dS=(2-fx7S I SHNIGSF=HACE—FxvS I] SHIGSIWIL/G=03939dS +vEOZ 
IdGOMI-VILSHL=VLISHL €C£O0Z 
7£02 S€E02 (IdDOML*° 39° VISHI)ISI 2EO2 
(AVLISO/XVLISOQI4NVIVtId=VL3SHL TEOZ 
ceO0¢d OL OD 
(AVLISO/XVLISG0)SNVLIVtIdGOML=VI3SHL OF0Z 
TEOZSOEO? (O0°O0°L9O°AVL1I3S0)4I €02 
G9O00TSEO? (AdSXVW°37° (SWIL/d)) GI 
(AG*AG+XG*¥Xd)4d1LHOS=0 202 
9O00TS 202 ( GdSXVW° 31° (3WIL/AQG) ) JI 
(AVL130) 4S8V=AQG 
AdW31-(E*S TSH) V=AVLI3—0 TO0Z 
9O0OT* TO? (GdSXVW°37° (3SWIL/Xd) D4] 
(XV¥V117130) 4S98V=XQ 
XdW31-(2*7*X) V=xXV1130 
LISOdN*T=X 9OOT OA: 
(CSCS T)V=HAdWSL SS (2° SIT) V=XdDW3IL 
LISOdN*T=I 9OOT OG 
, T+r=1 
VAN*T=F YOOT OA 


°OdSxXVW OL TWWNOZ YO NVHL SS37 SI d33dS JI 
(QASaer*s1)8 NI Y3S1N3 *T-LISYONST=F YOS T+r LIGYO 40 NOILISOd HDV3 
Ol fF LIGY¥O JO NOTLISOd HDV3 WONS G33dS GNV 3SYNOD *3ONVLISIG 3LNdWOD 


VUUYU 





606*976(VAN2OS°dN) JAI 
Q3IdNID0 SV3uvV JO YSEWAN SJNIWY3S130 OL 3SVHd LX3N OL LO3rOUd 


S46*906(0°0°19*% (7+ 1799 6G1)9) 41 
226*806 (0°0%19*%(T+1°9Erea1)@) Jl 
LISOdN*T=@I 226 OG 

O=NIGNSO=9VI4N 

LISOdN*T=9f 906 OG 

| y—dNxv=7 

; VINST=dN 906 OG 


°SNWNT0D AG SLNSW313 G3LVNOISSGNN 3AILISOd YO4 YDIHD 
te 109*+706 (0*%1L9°YS3YNNAN) 4I 
3NNILNOD 

(SN) MOYOINN 17VD 

VINS T=SN 266 OG 

O=YSSYNNN 

206*1S6 (ADVYLN°O3*VADVYLN) JI 

ae 3JNNILNOD 

106 *€06(0°03°3137390N) JI 

(dN) MOYOINN WIvD 

006*T06 (0°03°31L37350N) 4I 

(dN) IODOINN W1Vv>D 

VXIN*T=dN €06 OG 

YOVULN=VADVULN 


*SNOI131730 1LN3W3173 IVNOITLIGQV ON HLIM G3L31dWOD dOOT JNO TILNN 
SISVHd 11V 430 SMOY GNV SNWN10) NI SIN3W313 SNOINN YOsS ADIHD 


O=J31L13 130N 


O0=3 1D ADNSODHLOIMG=DHLGI MESOHLOIMG=HLGIMGSOH! d3GG=H1LdIGGSO=ADVULN 


°INNOD ADVAL JLVILINI 


S76 


806 


706 


cS6 


TS6 
£06 
T06 
006 


206 


UUY 


UUYU 


VUUYUYU 


UYU 


63 





SANILNOD ¢26 
JNNILNOD LO06 

SANILNOD O2L06 

€26 OL OD F T+NIGN=NIGN 0226 

CA0GlO econ Onaga OC UEen cal * ol) e—-(e+1*3i*Gl1)e8)3Sev)) 4] 00¢6 
DOC oCOecov\NEGiigagay (tt /+qea 1s] )G=-(C+7*31*81I8)3S8V))21 O26 
0L06*£0¢6 (H1d308°37°(((94+7S@1491)9-(2+75391°91)9)5S8V))41 616 
OZ0GcOLGUOTO 15" 01+ 1* ale ai en al 16 

LT6*OL06(GF°03°3T1)51 $ LISOdN*T=3I 0206 OO 


MOY LX3N OL SANNILNOD GNVY QG3SANO0S3SYNN SAVE OS JI 
°*IN3W373 Q3SLVNOIS30 3WVS SH!L OL INT LI3F Oud YIVd YSHLO ANV YOS ADSHD 


T=9OVISN 9T6 
EOG* OIG ple O°19°(8+1S9ISD1)8)51 F$ DI=QUNIGN $F GI=VUNIOGN 
T+NIGN=NIGN €T6 
ZUO- €16(DHLGIMG® 39° (0(L2+7°91SS1)8-(€+1S9rS8118)4S9V)) 451 O2TS 
OCT6SETEO(HLGIMG°37°(( (247991551) 9-(€+7°9F°G1)8)35S5S8V)) 41 2T6 
BOG Varo m EgaGg- 3 1°(((9+7*91*51)9—(2+1%9r* GI )8)4SGV))41 I16 
BOG* EUGr (0° 0 ree 06+ 3° dT * 1) 8) si 
LISOdN‘ST=DI L06 OQ 606 
eco, Ol OD S4NNITNOO 276 
SJANILNOD S26 
GL6*906(DHLOIME* SO? (0 (T-1S91S Er) 8—(e€+71*9r431)8)4SGV)) 41 626 
6L6*°906(HLOIME° 37° (1 (T-149T*9EF)G-(€+78Or*31)9)4S8V))41 826 
GL6*8L6 (H1d308° 37° (102-7991 °EF)G-(24+7°Er6391)98)4S9V) ) SI LLl6 
GL6SLLE6(O0° O° 19° (T+1S9FS31)98) 451 9L6 
9L6°SL6(GI°O3°3I)4S1 $F LISOdNST=3I1 G26 OO 
T=9V ISN 726 
Puce We. O- OL oo To] a) Gti 
JI=INIENSEI=YNIGNST+NIGN=NIQN 226 
MZoveocov snag 39 ((Ci—-1¢51*er)e—-(e+1*ar*gl)a)4aSev) 41 €L6 
CL6* CLE(HLOIMNG° 31° (0 (T-1S51S98F)8-(E+1°9Er*91)8)4SEV)) 41 6476 
176% 676(H1Ld308° a oe ate lod) G=tc+ 1°oh*agiid)3ISdv))4d] 8vé 
inc BevotOrO- 19° Ce— 1S oeor ) Gg) di 
LISOdN‘T=)1 L¥6 O00 946 


VUUUYU 





£26926 (1°O03*OW14N) AI 
3NNI LNOD 
ANNI LNOD 


Zuo Ooi nde 49° (4 (2+) °O1* >I] )g-(e+ 1*91*G@1 19) 43S9V) ) 41 
OSZ6S906(HLOIME°S7°( (024701491) d-(€+7S91601)8)5S8V) ) JI 
DuowmIcowmtl dag ad (9+) * qi «511 8—(c+1*91*01)S8)4S8V)) JI 

WuGeycou\ Oe O° LO° (le l*91°ah) aya 


, 
é 


OVT6SST6( 


LISOdN*‘T=91 47T6 OA 
€Z6SOvT6 (G1°O3°OI)GI 
0°0°19°(S+1*GI*51)9) JI 
LISOdN* T=dI OVT6 OG 
VYNIGN=GI ¢$ BYUNIGN=DI 
8€6 OL OD 
(dNSar*@l)LD39F0Y¥d W1VD 
JINNILNOD 

JNNILNOD 


ahe SOOCSHEGIMNd- 39° (€4+1*01 691) G-(T-15 51401 )98)4SaGV) ) SI 
6GT2*9O6(HIGIMG° 371° (((€4+7S0Q1S91)8-(T-1S91 40119) 4S9V) JI 
ST¢e*lte (HLd308°371°(0(2+7S01S91)98-( 2-19 916G118)45SaV) ) SI 


Gia s ict 


OF OO tea ale oi) ) aia 
LISOdNST=9I STZ OG 
Tineeal Cc (eh “U4 shat 


Auer erewo- O° LS* (e=1*SI1* diya) al 


°NWN 10D LX3N OL OD *3A10S34H LONNVD *0OS JI 


LISOdN‘ST=QI ¢T2 Od 
YNIGN=8I $ INIGN=1 
OG [ET Le so vAN-O4 > dN) aT 
906*OTZ(T*°O3°NIGN) JI 


*vauv agWwVS OL LD3rOXud 


GQINOD LN3W3145 Y3SHLO ANV GI 33S O1 LN3W313 LVHL ONINIVLNOD MOY 
JHL ADSHD N3HL SXD3HD NWN10D AG G3SIdND5DO0 GNNOS V3YV 3NO AINO JI 


= 


i 


Ov7T6 


716 


0S26 


S26 
726 


€26 
ST6 


OT6 


Oce 
cl2 
Ste 
61T¢ 
EVE 
Sid 


WIe 
etd 


i Oe 
Oi 


VUUY 





TE6*% C6 


826°6C6 


109*0S6(0°LD°YSAIYNNN) AT 
ANNILNOD 

(SN)JMOYOINN 1W1V5 

VANS T=SN £56 O*D 
O=YS3ZYNNIN 


Q1X3N OL 909 NOISSV $ SNNILNOD 
ye ebis: 


(WIVYLN®OS°VADVYLN) JI 
ANNI LNOD 

9¢€6%47€6 (0°03°313 1390N) JI 
(LN) MOYOINA 11V5 

Srot9eco, (0 OS" ala 13GN) a7 


(IN) JODOINA W1V) 

VANS THLN 726 OC 
ADVYLN=VAIVYILN 
(dN*GISOT) L539F0¥d 11V5 


°!>59r0Nd GNVY JLYNOIS3AG 


8€6 OL OD $F$ ANNILNOD 
AINNILNOD 
O°O==CII+ 1S Al* alia 
7*T=11 O€T6 O”A 


SCO Teo di ° Os ° Hie 41 

BO O° 25° C1 t 127) ot ) aps 
LISOdN*ST=AI TE6 OCG 
SNNILNOD 
JANTLNOD 
O°O-=(11+17S8FSA1I8 
7*T=11 0826 OG 

OEG*Ec6 Cal*OS“Al) 4a! 


(O°O0° Lo* (i+1*er sale) 31 
LISOdN*T=AI 826 OG 


N=(7+ 1989S @I)GSEC=(dNSN) HOVELS (84+ 1°EI SST) G=N 


*SINSW373 G3LYNOIS3G 


ao NWN10> GNV MOY WOYS SIN3W3739 YSHLO JAOwWSY GNV JLVNOTSISAC 


€S6 
S09 
906 
7&6 
9¢€6 
St&6 


B8e6 
L126 


Te6 
Of£ts 


66 


E26 
c£6 


8c6 
O8c6 


Of6 
6c6 


926 


VUUYU 








S996 05010" 07 19° (7+ 1* ar*e1) ey 41 
L=IVNgIN -L90 
69965296 -(0* O° HR eeaeele cee Orley edt 
T+NIGN=NIGN S96 
€66*S96 (DHLOGIMG® 39° (((L476a1691) -( 6475 Gwe essay) ) sie yoo 
Zoe Om encima tee’ Gl’ S1)e—-(e+ 1° er*dl)e)3Sevi) sl 996 
BoC Jou (Hi agiGd an (( 194 1°ol*>1)a—-(c+7°ar*al)ad)3sSev))4l €96 
eGo r9o0(G"’O°19°(S4+1*dI*>1)a) 4! 0296 
6296. 126610" 19* OV 1gN? a) 
LISOdN*ST=)I €66 OG 


°3SVHd LXJN NI ADNWVdNDDO VW3YV YO3S YD3HD 


O=NIGN 196 
Z66*196 (0°0°19°%(1T+17*9Er*agl1)g)4I1 
O=OVISN $ O=9OVI4EF 
LISOdN*‘T=GI1 266 OG S86 
O=OVI4NN $ O=NIQ7 
LISOdN*‘T=9f 066 OG 
T=9V14X ST+dNN=dWW $ 4—-dNNx=7 
| O*1T-=(Zx*Z7F*SZ1)D OTF 
€*T=ZyA Oly OG $O*T-=(7Fr*ZI1)A 
~elISGdN«l=z27 Oly Od $ LISOdINSI=Z1 OIL» od 
ONS T=dNN 066 OG 909 
~GLX3N OL O9 066 


°SAISNIDXS AVIVALANW LON 3YV SSYHL SOML *3NO S3ASVD 


*T+l JSVHd NI 
YIVd V OL SLOAFOYUd I ASVHd WOYNS NWNIOD DJI4ID3dS V WONS YIVd ON--YNOA 3SVD 
T+#f.3SVHd NI YIVd JYOW YO OML OL SLDIFONd I 3SVWHd WOYNS YIVd JNO-3IXH!L 3SVD 
°T+I 3SVHd NI YlVd 
VOL L53F0UYd I 3SVHd JO NWNI0D DJISID3dS V WOYNS YIVd JYOW YO OML--OML 3SVD 
Blt d5voo Nil divad 

NOWWOD V OL 1Dar0ud I 3SVHd JO 12V WOUS YIVd 3YOW YO OMI--3NO 3SVD - 
*SNOSV3Y WHIASS YOS LNIOd SIHL LV G3LVISYYNOINN 3G AVW SLISOd 


UUY 


VYVUUUYUUVYUUYUYUY 


_ 





1S1G=(@1*D1)a 
(14216 @r*@I)NIGLSIG=1SIG 
TL96*LL96 (0°0°19° (84791691) Q)4I 
(847° QI *D1)G=(ESEr*Ql ISDI=H(Z*régl DSLSIG=(T*Ersgl)d 
HL9O6STL96 ((TSEr*G1)d°399%1S1G) 41 

: T=9v14r 
7196*8196(0°0° 19° (2 EréGl DI SIS( TDI Ors GIINIGISIG=1S10 
OL96*TL96 (T*LO*NIGN) SI 


ig 


/ 


en, ee 


at eee cl aleve leech 
ONV IT+N LIGYO LV F dIHS N33ML39E NOTLVISYN woes G3LDSrOUd VSYV 4O YSLN3ID 
OL N LI@dO LV EGF dIHS JO NOILISOd WOYNS SDNVLSICG SSLVDIGNI (r*t)a 
°T+N SASVHd YOS YSEWNN dIHS SSLVDIONI & 
°*c+N SASVHd YOS YAGWNN dIHS S3LVDIANI I 
--NOTLVWYOSANI ONIMO11N04 SHL SNIVINOD (r*I)a 
"QSLVNOISSGQ AISNOIAZSYd AOVYL 4I NOITLVNOIS3IG S3LVDIGNI €=% 
*c+N S3SVHd YOS YSEWNN dIHS SAJLVDIGNI 2=% 
SlUNe edo shy bh dins 
GNV N LI@y¥O LV fF dIHS N33M13E NOILVISY WOYNS G3LDSCFOUd NIG 4O YSLNID 
, Ol ct+N 11990 LV (2@SFST)D dIHS JO NOILISOd WOYNS JDNVLSIG SSLVDIGNI T= 
°N SSVHd YOS YAEWNN dGIHS S3LVDIGNI & 
*T+N SSVHd YOS YSEWNN dIHS S3SLVDIGNI I. 
. —~-NOTLVWYOSNI ONIMONIOS SHL SNIVINOD (ASFEITID 
| _ “SLSIX3 ITI 3SVD °G AVY 3SN LSNW 3M *T LO NIOGN JI 
“ °15A NMONSA 


\ . 


9196 
eZ gS 


7196 
CL96 
8196 
OL96 
6996 


VUYUUYUYUUYUYUUUYUUYUYYUYY 


68 








; Lf 43LVNOIS3G 


*VaadV Gg3ld3arFO0Y¥d JO YSIN3SD OL SDNVISIG LSSTIVWS SVH 


O07 OL OOSIDI=ADIWSXC=HACWSAI=AIWSVIXAN OL O£66 NOISSY 
Xie Dot avi) ) ones TIVogec OX Al) D=DOIS$T=9OV 1458 


Cte xt NIE) 


SANTI LNOD 

MC=Xf S$ (TSME SAT) D=VOWSL 
Pecomoccouli Mir * ANI )D* a 1* VOWS Lad 

Ceoe seooul O=45° (Tt AC * ALD) 41 

LISOdNSXI=MF C66 OASAFHRXFS(IT SX ISAT DAVGWALSAC ERX) 
| C€66 OL OOSANNILNOD 

VeOouuZoe MUO TOPO CENT ONT) 2 al 

LISOdN* T=AF 7€66 OCGSLISOdNS T=AI O€66 OA 


j *LNS3W313 G3SLVN9ISSO V OL ONILDACOYNd SIN3W3139 3SOHL YO4 3AN0S3Y 


bs Cpls) 1). O 4 OoL 


G766°7266(T°39°9OV TANN) GI 
1266*066(0°03°NIQGT)SISSNNILNOD 
OVISAN+9OV TSANN=9V 1ISNN 
T+NIQGJ=NIG1 
(O° LO*NIYUN) JAI 
SANILNOD 
SANILNOD 
(8+7°Q@rI* OI )G=l(esais 91)9$91=(2¢415 GIIDISLSIG=H(T*3SI*QGI1)5D 
CTL6SCG66l(( TESTIS GI) D°S39°LSIG)GI 

T=9V1SF 

Qeeoenouo OI) (2°44) "oI )>ySsTst Teo) *S1* Sal YNIGISIG=LSIa 
766*TL6 (DHLGIME*S39°(( (247% @1S91)9-(€4+179391691)98)45SaV) 41 
OL6*TL6 (HLOIMG°37°( (0247591591) 9-(€4+17551 491098) 4S9V) AI 
766%696 (H1d3098°317°(( (947° G1IS91)9-(2+17S5I1*G1)9)45SSEV)) GI 
766*896 (0°0°19°(T+1S3S31°QG1)9) 41 
(@r°O3°SI)AISLISOdNST=31 4766 OG 


T66*TL66 


0896°766 


elt alo VHd wae 
LO4AFOYNd YIVd YAHLO ANV JI 33S OL I 3SVHd JO @I MOY YXD3HD 


euce 
€£66 
GS€66 


OL66 
7€66 


UU YU 


69 





me P_- 


00% Ol OOSVIXIN OL S766 NOISSVESAGXI=DJIWSEC=HACWSQXC=AIW 
(QXF*QXI Sar S dWWIDIS3SG W1V) 

007 O1 OD $ VIX3N OL S766 NOISSY $ DDI=ADIW $ DXF=ACW FS DXI=IW 
XE Ooo ox dws s0 aye 

(2S DXF SDXT)VD=991 

€766°¢C766 (QdW31°35 1° DdW3l) 4! 

O36° 2766 (OMNIS? Laan) a1 

JANI LNOD 

AF=QXFE AT=QXI $ CAL SAT) OG=QdW3l 

oa TH66(QdW31P 11° (AFSAIIQISAI 

t+ t340-) 448 

F Scoo U7G6(050° 39° Al °A1)G) 431 

LISOdNST=AI 6€66 O00 $ LISOdNST=AL 6€66 OG 

JANI LNOD 

ACF =DXF SAT=DXI $F (TSAPSATID=DdWISL 

Peoumucoo od Waleebnd= (T° ALS A1)D) 41 

it base as 

weco, Leoo "GO" Ondo CE * APY AL) ) J] 

LISOdN*T=AI 9€66 O00 $ LISOdN*‘T=AF 9€66 OQ 
O=LADISAdSXVW=AdWILSAddSXVW=IdW31L 


°1I 3S1VNOISSQ *3DNVLISIC 1S3M071 HLIM GQ OQNV D NI AYLNS YOS ADIHD 


JANT LNOD 

VIXSN OL O9 S$ O° T-=(IWSZF)Q 

O°T-=€2FSDIW)IG $ LISOdNST=Z2F 204 OQ 

 VWILX4dN O1 09 S O° T-=(2F*Z1)Q 

LISOdNST=ZI GO” O00 $ LISOdN*T=Z2F SO OG 

WO oecoyuedh Os (hates JniNTLNOD 

We eee Alm O° T=—(2°2@*71)) $ O° 1T—=(1*Zr*Z1)>D 
BOtvmeaO7s COIW Os ec* 26 * 21) 5) 41 

LISOdN‘T=Zr 904 00 $ LISOdN‘T=ZI 90% OdSO*T-=(Zx*SFWS ZF )D 
O*T-=(2AS ZFS IWIDSEST=aZy TO” OGSLISOdN‘T=7F TO OA 


*SNOTLV 13¥dOD 
Q31LVNOIS3Q 3SOH! YO4s G GNV D JO SNWNT0D GNV SMOY WONS SASTYLNS SAOWSY 


€766 


C766 


£766 


6€66 
T7766 


0766 
9€66 
8E66 


LE66 


S766 


UUYU 


O€66 
CO” 
70% 
SO” 
tO” 
90% 
L£0% 


TO” 
00% 


UUWY 





70 





~ 


“em. ene, 


JNNILNOD 

JNNI LNOD 

JINNILNOD 

-=(QNSON) XDVYULSI(GNSON) AXDVYL=(OGNSN) ADVYL $F LIGHYON*SGN=ON 80E OQ 
GOES LOE (I QNSN) XOVYLTOS® (GNSONINDVYL) SI 

90E*GOC(N*ORF° DN) ASI $ ADVYLNST=DN SOE OG 

~~ OOESHVOE(CLIGHYON*®LISGN)4I S$ JNNILNOD 

VN=QN 

ZOECSEOE(OP LOT (CYNON) XDVYULISI & LIGYON*ST=VN ZOE OA 
OOESTOE(OPLOS (TSN) HOVYLI AI S$ YOVYLNST=N OOE OA 
(9°64* =DHLOIMBHOT* XS*9°%64T 

‘6 =HLOIMQH6*XS*T°94* F=HILdSAQGH6SXS*ZI* H=JIWDADNHE6*XTT*//) LYVWYOS 
DHLOIMBSHLGIMGSHLd3S0G* 3 1TIDADNS T2476 INIUd 

T-31IDADN=S31DAON 

206 OL 09 $ HLOIMG-IdOML=DHLOIMG 

Id=HLOIMG 

1966S 0966(Id* LO*HLOIMG) SISHIGIMGxO*Z=HLGIMG $ H1Ld309x0*°Z=H1d3dg 
GO9 OL O09 $ GLXI3N OL 1566 NOISSY 

0276 *0SG66(2* 19°31 DADNI AI S&S T+31DADN=S 1DADN 


°1lV3d3uy¥ GNV 32ZI1S VSYV SSVSYONI 

Q3A10S3SY¥ LSA LON SYV SVSYV AldW3 OL ONILDSCFONd SNOILISOd SSOHL 
°VaeV AQSaLlDSrFOY¥d JO YSLN3D OL SADNVISIG ONIYVdWOD AG SLDITANOD V3uV 
INIAIOANI S3SVD SHL SA10SSY¥ OL GSADSHD N3SSG MON SAVH SXDZBHD LN3W3S13 
SNOINN AG QSA10SSY LON 3YSM LVHL YNOS—-3NO S3SVHd NI SLOVLINOD 4O W1V 


INNILNOD 
686*066(0°03°OV140) JI 


*ALdW3 3YV QO aNV D 


OO 
SOL 
SOL 
nO 
90t 
70 
2 10NS 
COG 
TOE 


Lea76 
L09 

02476 
T1966 
0966 
1TS66 
0S66 


066 
086 


VUUVUVUVVUYU 


UU 


Ge 





QNJ 

NYNL3aYy 

JANILNOD 

O=DJOINNN & (dWSErsGI)LI3FONd WIV) $ T=D0INNAN 
70S°O00G(0°0°19°l#+71°9OrSGIIGIAI $ N=GI 
OOSG*SECOSG(T*O3Z°W) SI $ ANNILNOD 

GI=N 

GOS*00S (T°1LO°W)SI $F Tt+W=W 

ROS = cos (0°0"19"*{14+71*9r*el)9) JI 

LISOdN§* T= gl TOS OG $ O=W ¢$ LISOdN*T=9fF 00S OG 
7-dW*7 = 1SO=JOINANGO=313130N 


DIOINOANS LISOdDNS GS DHLGIMET 


Sa Lalind * Bios ell YSAYNNN SOV ISAS VANS SLIISONS ADVYLSADVYLN NOWWOD 
: agesewAVII* (91° O1* 01 )9 NOTSN3IWIG 


°*IN3W3 13 
JATLISOd S3NO AINO SNIVLNOD HDIHM NWN10D ANV GNIS OL G SXADSHD 3JNILNOY SIHL 


(dW) IODDINN 3NILNOYUGNS 


ONS 

JANI LNOD 

((CXTSCII9EDSEXSCISXETS/) LYWHOS 

(LIGYOHS*XTT*/)LVWYOS 

(ADVYLNS T=NS (WEEN) XDVEL) SW SB8O0e LNIYd 
((C(XTSZ21)9€) ©XE* HOVYLHS*SXTT*//) LYWYHOS 

LIGYONST=wr 600€ OG $ LOOE INTYdS(ADVYLNS THAIS (I) L)*€00€ LNIYG 
T=(1)1 

SOVULNE T=I 900€ OQ 


SLUIASSY LNIYd 


pes. pyiv 


3%) 


opr 


00S 
70S 
€0S 
Tos 
SOS 
a (ONS 


O10€ 
2SI 

800€ 
LOOE 
600€ 
EOOE 
Z200€ 
900€ 
276 


UUUY 


UYU 


72 





GN3 

| ( (ONY) 4SODx(T 
9+ TS TSX)GH( 2419781) GKO%Z—-ZeH(9+TE TSH) GH+Zxe(Z479F51)G) 5SLYOS=NIGISIO 2 
DYV-IdOML=9uV T 

| 2*T(id°L9°9SUY) AI 

((L+791 69) G-( E4790 61)9) 5S8V=9UNV 

“GBTEBT°9=IdGOML $F LZ6STHI°E=Id 

DOINANS LISOdN* GS DHLGI MGT 

HIG IMG *H1GIMG*SH1Ldd0g* YSSYNNNS9OVISNS VINS SLITIGNS YOVYELSYDVYLN NOWWOD 

(S*SGZ2)XDVYLS(9TSOTS0T)G NOISNAWIO 


*HLIM Q3NY3DNOD 3YV 3M V3S"NV SHL 40 
YSLN39 SHL OL NOITLISOd G3I41ID3dS VY WOYNS JDNVLSIG S3LNdWOD NOILONAS SIHL 


(TSA STINIGLISTIO NOILOINAS 


ONS 
Nanlsay 
O=DOINAN S$ SJNNILNOD O08 
(dWSar’@l)LD35FO0Y¥d WWD $ T=D0INAN 4708 
708*O008(0°0°19°(7+18Er*GI)G)41 $ N=Or €O8 
O0O8*eO8 (T°O03°W)ISI S$ ANNILNOD TO8 
8r=N SO8 
GO8*O008 (T*LO°W) SI S$ TH+YSSYNAN=YSSYNAN 908 
Pee Ue CO SO TO ttt er *adlyeisl $ LtWeN 2og 
rOe coe CO°O°1OS°(1+1*er*gI)e)ai 
LISOdN*T=E8r TO8 OO $ O=W §$ LISOdN*T=8I O08 OG 
7-dW* 7= 1SO=SOINANSO=3L135 150N 
DOINANS LISOdNS GS DHLGIMET 
*HLOIMGSH1d308* YS3SYNNN SOV ISAS VANS SLITSON SAXDVYLSADVYLN NOWWOD 
(S*°SZIADVALSI(OTSOLTSOLIG NOISN3WIG 


°1NaW373 
JAILISOd 3NO ATNO SNIVLNOD HDIHM MOM ANY GNI4 OL @ SYD3HD ANILNOY SIM 


(dW) MOYOINN SNILNOYBNS 


VUUUWY 


ve, 


Wig UU 








SNNILNOD 

dI=INIGNS DI=YUNIOGN 

T+NIGN=NIGN 

009* E09 (DHIGIMG* SON (( (L496 gI S91) G-(E+1% drs Gel araody ) 741 
UO cOgtniGina Ss yeii (itn dgl* ol) g=te+1% ar *d})9)4S9v))41 
009* ZO9(HLd390g°31°(((947SdIS51)9-( 24+ 18 df*dI)G)4S9V) DAI 
a T+W3 1S3N=W3 14N 

OOF OeICO"OS1S*14+ 1)*gl<51) 4) a1 

LISOdN*T=)I 009 OQ 

079*S09(7°03°%ddW)4dI 


°3SVHd LX3N NI V3YHV G3IdNDDO YOS ADZBHD 


O=W313N 
7-VdWe7=1 6 O=NIUN 


°AIVEG SSLVIIONI T=9OVISA *GevVMYOS NOILDIBCONd SALVDIGNI O=9V144 


O=9V 1SAASddW=VdWSEF =df SAI =d! 


DJOINANS LISOdNS GS DHLOIMGT 


*HLGIMG*SHLdIjdg* YSAYNAN SOV 1SNS VANS SLITSONS ADVULSAOVYLN NOWWOD 
(S°SGZ)ADVYLS(9OTSOTSOTIG NOISN3WIG 


°3ISVHd LX3N NI G3IdNDDO SI V3YNV SI 3NIWYN3L3O OL SLYOd|aY NOILISOd 
JO SYIVd JLVGIGNVD G3I4IDadS WONS YDVE YO GV3HV S1D3CONd 3NILNOY SIH 


(ddW*@P*@1)L>39C00YNd ANILNOUBNS 


009 


E10), 
0209 
COF 
t09 
Oc9 


O79 


UUYU 


A, 


U 


UYUY 


74 








ElorcO lt Oo 

INIGN=df S$ YNIGN=dI 
CT9*SSO9 (VAN°OSR°VdW) 4s! 
I+VdW=VdW 


Sslovidees Vile O se Loa odd 


L=D0INNAN 
(DISdI*dffvdW)9ISAG W1VvD 
SO9*TI9 (Z2*O3*NIGN) SI 
ANNI LNOD 

ANNI LNOD 

L+NIGN=NIQN 


ae Ole WORM Aye Crea edi’ Sl} e-(e+ttai*dal )g)4Sev)) 41 
OOO wre i Mesa. (02+)? dl Seiler 3° 31°O1)a)489V) 141 
2Og Gog vhi added sn (((94+71*°d1* 51 )e-( 2+ 1831 SGI 19) 4S8V) 14! 


Boge O90" O° lO" (lL. )}*a1+ dl) dad) al 
LISOdN* T=]3I 909 OQ 

O902. 209(0° O° 1° 16+ 1*Gh*)1)9) ai 
LISOdN* T=GI 0909 OQ 

YNIGN=)] 

TT9 OL O9 $S$dI=INIGNSOI=YNIGN 
ANNI LNOD 

Wel oco( Ue O19" (St 1sdl*aI)g) 41 
piSOdN> b= 1° 769-00 

SOC Cees tees DOINiINI al 
GO09°2eo Ghigo Wadanoo 

eee ZO gee) 4° NEN ea) 


Paw NOTSaG ON Sie °NOILISOd SWVS OL LO3SrOs8d NV) 


MOVYL YAHLONV JI ASS OL ADZSHD NAHL *VANV NI NOILISOd ASANO AINO JI 


Cg 


2S9 


Lt9 


0909 
909 
Oe 
0609 
Ole) 
809 


£09 


709 
OG 
oe 7 


Boo 
Crs 


beg 


UUYU 


UUY 


75 








74> 


ON3 

' N&UNLSY $F$ O=DJOINNAN 

829 OL OD FO=WATANSOINIGN=EdFSUNIOGN=dI 

ECO STOW Os" vaw) sl 

T-VdW=VdWS T=D0INNN S(DNIOENSdI SdfSVdW)OIS30 WIVD 
ST9*9¢9(2Z°OS°NIQGNI SI $ ANNILNOD 

SNNILNOD 


T+NIGN=NIQN 


Te9 *S29(DHLOIMG* 39° (( (€+7°01S31)9-(1-71591 401 )98)4S8V) DAI 
OUNCO*SCO(HIGIMG° 371° (( (E+ 7SO0I S31) Q-(T-1S91*01)9)45S9V))41 
Te9* 9e9(Hld309°37°(((2+7501531)8-(2-15914G1)9)45SaV)) 41 
reo secs § 00° 0" 19° (1+ 1* Ors a1 ee 3] 

LISOdNST=31 [29 OQ 

Oe cc O° 0° 9*4e-1*51*G1)qg) 51 


LISOdN*‘ T=GI OT29 OQ 
INIGN=)!1 

929 OL OD SDI=DINIGN F dF=YNIGN 

SANT LNOD 


eee t7ot OP OeTo c= 1* D1 * df )a)si 

LISOdNST=)1 6€9 OQ 

SGT9*S8E9(T*POS°0INNN) SI 

ST9*LE9(T®O3°W313N) JI 

9E9SOCI(T*O3S°NIGNISI S$ SNNILNOD 

JI=OINIENSdF=4NIGN S$ T+NIGN=NIGN 
ITI*®6TI(DHLOIME* 39° (0 (T-9591% dr) G—-(€+ 1% dr édl)a)45S9V) )4I 
O61T9*6TO(HLOIMG 237° (((T-19916 dF) a-(E€+ 1% dF %dI1)G)4SaV) )4I 
eee ed dae alettico1* D1 *de)d—-(c+ 1*dr*d1)9g)4Sav)) SI 
T+W3 1T3N=W3 713N 

Pie ac ot US Oso: Ve 1 Oo) * df )g).51] 

LISOdN*ST=5I1 9T9 OG 

9-VdWx7=1E O=NIQN 

O=W3 ISN $ T=9OVISA S$ Af=adf $F$ AQl=dI $ ddW=VdW 

SIS° 91S (1 * io ddW) SI 


*SSSVHd Y3SIT8V3 OL ADVE LD3FONd *3NO NWHL Y3LV3YNO 3SVHd JI 


ST9 
ECS 


9729 
Glce 
boo 
SZo 
0729 
idee) 
Ea 


Cag 


Oc? 


c79 
Ged 


76 


8E9 
LE9 
969 : 
919 | 
619 
0619 
8T9 
LT9 


“GE9 


89 
719 
S09 


VU YU; 








N=(7+7°9F*SQGI)Q 

GF=a(VdWSNIADVAL 

GI=(T+VdWSN)ADVYL 
DI=a(2+VdWSN)ADVAL 

N= ore a oe i) a 

AIVYLN=N 

T+AVYLN=AADVYLN 

| Ole 2 Lao. 
: T+A=A 
(8+ 1°SI*S1)G=N 

SOE SOL COO 1D (s+ 1*uls ot) a)al 

O¢L OL OD SEF=ACFII=AIST+VdIW=AVdGWST+AR=n 
POTESSOUVO* OS * 3) Sl 

Ria ea to 17.3 

DT =a(Ct+VdWSN) AVAL 

GI=(T+VdWSN) AVAL 
O-=(c+VdWS NN) XDVUL 
O-=(T+VdWSNN)XDVUL 

(8+ 71°@ISOD1)G@=NN 

VO eos" OL" (8+ 1°8l*d1)a)al 

(7+ 1°9@F SQI)G=N 

COMO OOS 17+ 1°8r*e17d) a1 


OSLSOOL(0° 03° OV 14) SIS 7-VdWev= 11SEC =F SEI =1S0=ASGdW=VdW 
SOINANS LISOdNS 8S DSHIGIMGT 
*HLOIMESH1d308* YSSYUNNN SOV 1450S VANS SLITSONS ADVULSADVYULN NOWWOD 


Woven dso aN T ead el 
CIS cum oVvaol (91 SOT OL). NOISNIWNIG 


*SIN3W313 GSLVNOIS3SQ SHL 3O NWNT0D ONV MOY SHI WOYS 
SINSW313 YSHLO S3135173530 GNV SNOILISOd G3I41IDSdS SSLVNOISSG 3NILNOY SIHL 


(DIS@I*S ars ddW)9OIS3G0 JNILNONENS 


OTL 


60L 


BOL 
COL 
GOL 
SOL 


VOL 


col 


Los 
OOL 


UYUYU 








GN3 
| Nunisa LOL 
O€L OL OOSEI=A(T+VdWSN) XHOVYLSACHAIVdWSN) ADVUL 
N=(7+7°9F SQI)GSN=(799I1*S9FIG $ DI=H(T—-VdWSN)IADIVYL LSL 
ADVYLN=EN $ T+XDVYLN=XDVYLN 962 
LGL OL 09 $ TtHX=0 $ (1SDI1SEC)G=EN SSL 
9SGLSGGL(0°0°9LO*% (1991S EGrIG)4I 2S2L 
O€L OL OOSDI=APSEL=AIST-VdW=VdWSTt+X=AD . 
N= (9891*°98F)98 $ DI=HCT-VdWSN)XDVYL VSL 
LOL OL 09 $ N=(1SDI1*9FIG $$ GF=A(VdWSNINDVYUL $F$ DI=A(T-VdWSN)XDVYL 
O-=(VdWSNN) XOVYUL S$ O-=(T-VdWSNN)IXDVYL SF (TSDI*SEC)G=NN ESL 
HGLSESLIOPO° LOZ (TSDISGFIGI4AI & (7+1°EFSGI)G=N TSL 
/ CGLETSGL( O° 0° 19° (H+ 7°EGrSGl1)G)4I OSL 
O€L OL OOSDI=APSECL=AIST-VdW=VdWST+X=x TSOL 
LOL*UGOL(0°035°X) 4I OGOL 
0S602*S02(0°03°9V144) 41 
SNNILNOD SEL 
3NNILNOD 6€EL 
0°0-=(NITSF SOIL) ASItHI=ANT 
7*ST=HI 6€Ll OA 
T=3L37350N BEL 
GELSBEL(OPO*LOT(TH7SC SOL) GSI LEL 
LELSSEL(1°O3°OI1)GI 
LISOdN*T=9I SEL Od 
3NNILNOD TEL 
3NNILNOD vel 
0°0-=(AX1S9OI*S 1) QS 04HI=N7 
ST=HI vel OG | 
T=3L13730N €EL 
eelLsTel (fF °O3°OI) SI ZEL 
TELPZEL(O°O°LO®° (THIS 9I1 $1) G)41 
LISOdN*T=9I Tel OA 
y-VdWxv=1 O€L 


78 








00€00000 
06Z00000 
08Z00000 
01200000 
09Z00000 
0S¢Z00000 
0%Z00000 
0€200000 
0ZZ00000 
OTZO00000 
00200000 
06T00000 
O8TO00000 
OLTOO0000 
09TO00000 
0STOO000 
0OvTOO0000 
O0€TO00000 
OZTO00000 
OTTOO000 
0OTOO000 
06000000 
08000000 
0L000000 
09000000 
05000000 
07000000 
0€000000 
0z000000 
01000000 
00000000 


JivVIAZG 3AILVO3N YOS ADIHD 


v¥ Lvol4 
93N 31 LNSW31dWOD 


Mer elels 
H194% BSHL OL 2 LOVYULENS 


94W INdNI JO SSIYCGGV=*x 
Guev 3HL OL 2dOWw 39NGsY 


9Y4v LAdDNI 3O SSSYCQV=Hx* 


texas 24s 


SINAN OUV LSYI4S Las 
JINAN OV LSYI4d LAS 


ao ay SdNO2 J oa0 dis 


XSQONI SAWS 


€+SNOIAQY 
3SVuy3 

9¢ 

@) 

Ly 

Q@l1S02 
3SVuy3 

I+x 

3SVUY3 
SNODAQY 
JIAGY 
3SVy3 

HE 
T+SNODAQGY 
ETSAIgS 

a 

ST 

3SvVay3s 


Pixs 

T 
T+d1AQuy 
d1AQGa 
72 
t=L1x3s 
0 


E9ASGNY 


L1X3 

%% 
EQAZOGNY 
E9AZQNY 


ce 


WIS 
Aso 
S711 
VN3 
Su 
VOS 
IV 1 
dv 
WLS 
uMs 
df l 
Qva 
WLS 
AOS 
TNW 
VQ) 
INS 
Ess 
YN3 
ae es 
IN] 
WS 
AWS 
Suv 
Ws 
VO 
A 


ATS. 


‘om 


AMINA 
LNIQI 


dIAQY 


CIAZINY. 


tie, 





02700000 
0T7O00000 
00700000 
06£00000 
o8e00000 
OLE00000 
09€00000 
0S€00000 
O07€00000 
0€€ 00000 
0200000 
OTEO00000 


yY ONV GNODISS JO SSAYOCV = x# 
NOILVIA3G GYVGNVLS Ag JCIAIG 


O- 
aie Vee Lo 7 Cui 
OODDDD000000009L 
O000000000000002 

GZtTcOLlodecl 
t 

HF 

% 


XH 


C+SNODAGY 
€+SNODAGA 


SINI4 
ONJ 
ON3 
L150 
JEDI. 
Ze 
LO 
IAG 
Sod 
PUSS 
INA 


eee 


AW 4 
WIS 


mo ox a 


SNODAGY 
ela le 
aS Va 


LIX3 


80 





APPENDIX II 
INPUT AND OUTPUT ~ DETERMINISTIC MODEL 
The variable names used by this program are defined in the program 
listing contained in Appendix I. No provision is made for reading: the 
values of parameters into the program. Changes in parameter values are 
made by inserting punched cards into the program deck prior to compila- 
tion. The parameters which determine the formulation of the problem 
are; 
1) NPOSIT 
2) NORBIT 
3) SIZESQ 
4) CUSDEV 
5) SPDDEV 
6) MAXSPD 
7) MINSPD 
8) MX 
9) NUNIF 
10) NPROB 
11) BDEPTHO 
12) BWIDTHO 
13) BWIDTHCO 
14) TIME 
The DIMENSION statement for all arrays must be changed if the value of 
NPOSIT exceeds 10 or if the value of NORBIT exceeds 5. The values of the 


other parameters may be changed freely. 





The sample output which follows illustrates several points: 

1) Governing parameters are printed at the beginning of the out- 
put. The number of course changes and the number of speed changes to 
be made for orbit J, as determined by random numbers, are indicated. A 
course change or speed change indicated for orbit J determines the course 
or speed for the period between orbit J-1l and orbit J, 

2) The column labeled TRACK gives the position of each ship at the 
time of the orbit number listed in column TIME. The course and speed 
listed is the course and speed for the ship indicated for the period 
prior to the time of the position report, 

3) The course, speed and position for each ship at the initial 
orbit is generated by use of random numbers. The program then determines 
how many changes of course and speed will be made, which ship will make 
them, and the sign and magnitude of the change. The later courses, speeds, 
and positions are then computed and used as a basis for correlation. 

4) NCYCLE indicates the number of times the dimensions of the 
constructed areas were increased before the final correlation was made. 
Final dimensions are indicated. 

5) The print-out of tracks indicates which ship was designated as 
a specified track for each orbit. Any track which contains the same 
ship number for all orbits is considered to be perfectly correlated. In 
this example Track 2 and Track 14 were erroneously projected back to 
orbit 1. The cause is not apparent without plotting all tracks. 

6) Note that a "mistake" in correlation as indicated in this 
example does not necessarily imply faulty correlation logic. Many cases 
of this type result from crossing tracks and erratic ship maneuvering. 


The "correct" logic may pick the "wrong" ship. 


82 


ee I eee 
SS yy» = 
bic © (a ae 


- a ee oo 


eee ee) 








7) The tracks containing all zeroes result from duplication of 
designations. For example, ship number seven may have been designated 
as Track 3 for orbits 3, 4, and 5. It may have been later designated 
as Track 6 for orbits 1, 2, and 3. The duplication is resolved prior 


to printing the results of the correlation. 
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APPENDIX III 
FORTRAN 63 PROGRAM = PROBABILISTIC MODEL 
The input to TRACK IV is a magnetic tape prepared by some program 

such as READOUT. This input tape may contain simulated position re- 
ports for one or more readouts. The data from one readout can contain 
from one to 100 position reports and must be written with a FORTRAN 
'WRITE TAPE i,L' statement where L is an array dimensioned as (100,3). 
Each position report occupies three words of the array L as follows: 

L(i,l) = Time of report 


L(i,2) = Latitude of the position 


i 


L(i,3) = Longitude of the position 

The array L must be TYPE REAL, that is all data must be in floating 
point. The data must be in the units of measurement described in the 
next paragraph. A problem termination is indicated on the input tape 
by an END OF FILE. No END OF FILE should appear between readouts. 

The method and units of measurement for the quantities processed 
by TRACK IV are as follows: 

1) Time — Time is measured in hours from the beginning of the 
problem to the end. Time need not begin with 0.0 for the first report. 
In fact, if the program READOUT is used a time less than 2.0 will not 
be processed. Time is a floating point quantity. 

2) Latitude - Latitude is measured in radians from 0.0 radians 
at the north pole to T{ radians at the south pole. It is a floating 
point quantity. 

3) Longitude ~ Longitude is measured in radians from 0.0 radians 
at the Greenwich meridian eastward through 27T radians. It is a float 
ing point quantity. 


4) Course = Course is measured in radians from 0.0 at north to 
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the right through 27{ radians, It is a floating point value. 

5) Speed - Speed is measured in radians per hour. It may be 
converted to knots by multiplying the value by 3436.6. Speed is a 
floating point quantity. 

6) Distance = Distance is measured in radians and may be con- 
verted to nautical miles by the same conversion factor specified for 
converting speed to knots, It is a floating point quantity. 

7) Probability - Probability is handled in two forms. In most 
computations it is a floating point value between zero and one. It 
is stored in lists by multiplying the value by 1000 and converting 
to integer form. This permits the quantity to occupy 12 bits in 
memory. Some computations in TRACK IV handle probabilities in their 
integer form. 

TRACK IV utilizes a permanent record tape on which to store track 
histories between readouts. The tape also contains the next available 
track number and a list of any uncorrelated position reports. The 
format of these lists is as described below for the arrays TRACK and 
OLD DATA. This tape is set to an initial blank form by TRACK IV before 
processing the first input from the input data tape. If it is desired 
to continue a problem on different runs, that part of the program which 
initializes the permanent memory tape must be removed before the second 
rule 

Two tape units are used by TRACK IV as temporary storage. These 
are identified as units two and three in the program. The input tape 
is designated unit one and the permanent memory tape as unit four. 

Data design for TRACK IV is as follows: 


1) OLD DATA and NEW DATA are arrays dimensioned (500,3) and 
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(100,3) respectively. Both contain position reports in the form: 
NEW DATA (1,1) = time of report 
NEW DATA (1,2) = latitude 
NEW DATA (1,3) = longitude 
2) TR LIST and TR VEL are equivalent names for the same array 
which is dimensioned as (31,4). Data is addressed as TR LIST if it 
is an integer value and as TR VELif it is a floating point quantity. 

a) TR LIST(1,1) bits 0 - 5 contain the tentative part of 
the track number. Bits 6 — 23 contain the basic track number. Bits 
24 ~ 47 are unused. 

b) TR LIST(1,2) bits 0 - 11 contain the probability of 
termination computed on the last sweep. Bits 12 = 23 contain the 
track reliability assigned to the track on the last sweep. 

c) TR VEL(1,3) and TR VEL(1,4) contain the time of initial 
track contact and the time of last extension respectively. 

d) TR VEL(2,1) and TR VEL(2,2) contain the latitude and 
longitude of the last known position respectively. 

e) TR LIST(2,3) is set to 777 octal if the list concerns a 
previously uncorrelated position report instead of a previous track. 
This flag indicates that there is no history for the track. 

f) TR VEL(2,3) and TR VEL(2,4) contain the last known 
course and speed respectively if the list concerns a previously known 
track. 

3) The remainder of the TR VEL/TR LIST array contains possible 
correlation in the form: 

a) TR LIST(I,1) = new data list index corresponding to the 


location of the position report in NEW DATA. 
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Die in Licht 2) Pg OF Po for the correlation listed 
until the list is normalized; MAPS 41 4? thereafter. 

c) TR VEL({I,3) and TR VEL(I,4) contain course and speed 
respectively as indicated if this correlation is a valid one. 

4) PL LIST is dimensioned (100,29) and is TYPE INTEGER; that is, 
data in the list is integer rather than floating point. In this array 
all correlations which involve a specific new position report are 
tabulated. Each entry requires one computer word. Bits 0 ~ 11 contain 
the Pp, or P59 46 appropriate untii the list is normalized; Lai yA) 
thereafter, Bits i2 = 17 contain the track tentative number and bits 
18 = 32 the basic track number for the track with which the position 
report may be correlated. Bits 33 = 47 are unused. An entry addressed 
as PL LIST(I,N) refers to the Ne possible correlation for the ite entry 
in NEW DATA. 

5) TRACK NR and TRACK are equivalent names for the same list. 

The list is eight computer words long. TRACK contains the permanent 
record form of the track history as it is stored on the permanent 
memory tape. If a data word in the list is an integer it is addressed 
as TRACK NR; if it is a floating point value it is addressed as TRACK. 

&) TRACK NR(1) contains the track number and tentative 
number in bits 0 = 20. Bits 21 - 47 are not used. 

b) TRACK NR(2) bits C ~ 11 contair the probability of termin- 
ation computed at the time of the last extension. Bits 12 - 23 contain 
the track reliability. Bits 24 = 47 are not used. 

c) TRACK(3) to TRACK(8) contain, in the order listed, time 
of track initiation, time of last position, latitude of last known 


position, longitude of last known position, last known course, and last 





known Speed. 

Table 14 is a sample output from TRACK IV. The numeral to the 
left of the radix point in the track number column is the basic track 
number. Ideally it would remain unchanged from the time a track is 
initiated until the time it is terminated. The digits to the right 
of the radix point are branch identifiers. They are assigned at each 
readout to identify all possible correlations for the basic track. 

Arrows to the left are not part of the print out. They were 
inserted by the author to show which of the possible extensions is 
the true one. The sample shown was the output for the test described 


in figure 12 for the four hour interval run. 
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APPENDIX IV 
FLOW CHARTS FOR THE PROBABILISTIC MODEL 

The flow charts in the following pages are intended to give the 
reader a general picture of the manner in which the logic of the 
probabilistic model was implemented in TRACK IV. Detailed flow egharts 
may be obtained by contacting the Operations Analysis Department of 
the U. S. Naval Postgraduate School. 

The quantity L,,, which appears in the flow charts is defined 
as the probability that a tentative track A.n established on sweep i 
is the true extension of track A in view of the possible extensions 
determined on sweep itl. I,,, iS computed as the maximum PsQ over 
all possible extensions of track A.n multiplied by the track reliabil- 
ity of track A.n computed on Sweep i. 

Lee is used to eliminate tentative extensions established on 
Sweep i which are duplicates of tentative extensions for other tracks 
established on sweep i. The duplicate with the highest Lax is re- 
tained. U,,, iS also used to determine which branch of a split track 
should retain the identity of the track and which should be redesig- 
nated. The branch with the highest LL, retains the former designation 


and all others are redesignated. 
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